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ABSTRACT 


This  thesis  provides  a  comparison  of  the  accuracies  of  two  methods  for  fitting  a 
WeibuU  distribution  to  a  set  of  aircraft  engines  time-between-failure  data. 

One  method  used  is  the  Maximum  Likelihood  Method  and  assumes  that  these  engine 
failure  times  are  independent.  The  other  method  is  a  Modified  Method  of  Moments  procedure 
and  uses  the  fact  that  if  time  to  failure  T  has  a  WeibuU  distribution  with  scale  parameter  X  and 
shape  parameter  P,  then  T*^  has  an  exponential  distribution  with  scale  parameter  The  latter 
method  makes  no  assumption  about  independent  failure  times. 

A  comparison  is  made  from  times  that  are  randomly  generated  with  a  program.  The 
program  generates  times  in  a  manner  that  resembles  the  way  in  which  engine  failures  occur  in 
the  real  world  for  an  engine  with  three  subsystems.  These  generated  operating  times  between 
failures  for  the  same  engine  are  not  statisticaUy  independent.  This  comparison  was  extended  to 
real  data. 

Although  the  two  methods  gave  good  fits,  the  Maximum  Likelihood  Method  produced 
a  better  fit  than  the  Modified  Method  of  Moments.  Explanations  for  this  fact  are  analyzed  and 
presented  in  the  conclusions. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  aU  cases  of  interest.  While  effort  has  been  made,  within  the  time 
available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors,  they  cannot 
be  considered  validated.  Any  application  of  these  programs  without  additional  verification  is  at 
the  risk  of  the  user. 
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EXECUTIVE  SUMMARY 


This  thesis  provides  a  comparison  of  the  accuracies  of  two  methods  for  fitting  a 
WeibuU  Distribution  to  a  set  of  aircraft  engines  time-between-failure  data. 

One  method  used  is  the  Maximum  Likelihood  Method  and  assumes  that  these  engine 
failure  times  are  independent.  The  other  method  is  a  Modified  Method  of  Moments  procedure 
and  uses  the  fact  that  if  time  to  failure  T  has  a  WeibuU  Distribution  with  scale  parameter  X  and 
shape  parameter  P,  then  T*^  has  an  exponential  distribution  with  scale  parameter  The  latter 
method  makes  no  assumption  about  independent  failure  times. 

In  a  fleet  of  fielded  aircraft  engines,  engines  are  repaired  when  a  fault  is  detected.  In 
addition,  some  critical  components  are  inspected/repaired  after  a  designated  number  of  flight 
hours.  When  faults  are  detected  the  engine  is  repaired  and  the  component  is  replaced  on  the 
engine  unless  the  engine  has  accumulated  a  total  number  of  flight  hours  near  to  the  hard 
inspection  time  of  a  critical  component(s).  If  this  is  the  case  the  critical  components(s)  are  also 
inspected/repaired.  In  nearly  aU  cases,  the  engines  are  not  completely  overhauled  to  restore 
them  to  new  condition.  After  repair,  components  in  the  engine  have  mked  ages  with  respect  to 
flight  hours  with  some  components  new  or  repaired.  Consequently  the  times  between  failures 
on  a  specific  engine  are  not  independent. 

A  comparison  is  made  for  times  that  are  randomly  generated  from  a  program.  The 
program  generates  times  in  a  manner  that  resembles  the  way  in  which  engine  failures  occur  in 
the  real  world  for  an  engine  with  three  subsystems.  These  operating  times  between  failures  for 
the  same  engine  are  generated  in  a  way  that  they  are  not  statisticaUy  independent.  This  is 
accomplished  by  using  a  program  in  S-Plus  which  simulates  failure  times  for  an  engine  with 
three  subsystems.  The  failure  time  distributions  of  each  one  of  these  subsystems  are  WeibuU 
with  fixed  parameters  determined  a-priori  and  based  on  real  world  experience.  The  simulation 
gathers  three  failure  times  for  each  engine.  The  study  includes  failures  for  two  different  sets  of 
20  and  100  engines  respectively.  Once  this  comparison  was  completed,  the  study  was  extended 


to  real  data  using  TF-34  engine  operating  times  between  failures  from  the  Naval  Aviation 
Logistics  Data  Analysis  (NALDA)  data  base. 

The  estimation  of  the  parameters  for  the  WeibuU  Distribution  which  fits  the  failure 
times  was  done  using  spread-sheet  software  programs.  The  accuracy  of  each  method  was 
measured  using  the  Kolmogorov-Smimov  test  statistic,  embedded  in  the  S-Plus  software 
package. 

Although  the  two  methods  gave  good  fits,  the  Maximum  Likelihood  Method  produced 
a  better  fit  more  often  than  the  Modified  Method  of  Moments.  Explanations  for  this  fact  are 
analyzed  and  presented  in  the  conclusions. 
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1.  INTRODUCTION 


A.  BACKGROUND 

1.  Engine  Reliability  Analysis  Program  (TRAP) 

The  Naval  Air  Warfare  Center  Aircraft  Division  Propulsion  and  Power  Engineering 
(NAWCADPPE)  was  tasked  with  analyzing  aircraft  engine  failure  data  to  estimate  aircraft 
engine  reliability  and  maintainability  factors  and  to  assess  the  effectiveness  of  existing  naval 
aircraft  engine  maintenance  policies.  The  ERAP  program  was  established  to  help  accomplish 
this  task.  It  generates  nine  reports  designed  to  identify  problem  areas  and  trends  in  engine 
logistic  support,  time  between  engine  failures,  engine  availability,  engine  readiness  degraders, 
and  other  engine  performance  factors.  The  data  base  used  in  the  ERAP  program  is  the  Naval 
Aviation  Logistics  Data  Analysis  (NALDA)  data  base.  These  reports  provide  summary 
measures  such  as  mean  flight  hours  between  engine  failure  and  other  averages  and  trend  lines 
on  these  averages.  They  cannot  be  used,  for  example,  to  compute  residual  reliabilities  of 
specified  additional  lifetimes  for  engines  with  known  number  of  repairs  and  operating  times 
since  last  failure.  Nor  can  they  be  used  to  determine  optimal  maintenance/replacement  policies 
for  critical  engine  components.  A  more  detailed  analysis  of  the  engine  failure  data  in  the 
NALDA  data  base  is  required  to  make  these  types  of  assessments.  Michael  R.  Caudill  [Ref  1] 
and  John  A,  Malsbury  [Ref  2]  made  more  detailed  analyses  of  TF-34  engine  performance  in 
their  Master  of  Science  theses  than  those  provided  under  the  ERAP  program.  Mark  E.  Mlikan 
[Ref  3]  developed  an  integral  equation  for  determining  optimal  no-build  times  for  critical 
components  of  TF-34  engines  relative  to  their  hard  inspection  times.  All  three  of  these  theses 
used  a  Weibull  distribution  to  model  the  probability  distribution  of  the  flight  hours  between 
engine  failures  for  the  TF-34  aircraft  engines  archived  in  the  NALDA  data  base. 
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2.  NALDA  Data  Base 

NALDA  is  an  operational  automated  information  system.  It  is  a  centralized  Naval 
Aviation  3M  data  base  supported  by  the  Defense  Megacenter,  Mechanicsburg,  Pa.  A  user  of 
the  NALDA  data  base  must  make  coded  queries  to  the  system  to  obtain  specific  information. 
This  requires  considerable  experience  with  NALDA.  NALDA  contains  numerous  data 
elements  on  Navy  aircraft  engines  by  engine  type  and  serial  number.  Anytime  an  aircraft 
engine  is  inspected,  maintained,  removed  from  an  aircraft,  repaired,  transported,  reinstalled  or 
discarded,  a  report  is  submitted  that  is  translated  into  the  NALDA  data  base.  The  NALDA 
data  base  must  be  “cleansed”  before  it  can  be  used  to  perform  some  detailed  statistical 
analyses  such  as  fitting  a  probability  distribution  to  the  empirical  cumulative  probability 
distribution  of  operating  times  between  engine  failures.  M.  Caudill  did  this  cleansing  of  the 
NALDA  data  for  TF-34  engines  provided  to  him  by  support  staff  at  NAWCADPPE.  Details 
of  this  cleansing  process  and  its  purpose  are  described  m  Caudill’s  thesis  [Ref  1]. 

3.  Operating  Times  Between  Engine  Failures 

The  NALDA  data  base  reports  the  operating  time  between  engine  removals.  Engines 
can  be  removed  for  non-failure  reasons.  Caudill  [Ref  1]  established  a  data  base  of  operating 
time  between  engine  failures  for  TF-34  engines  by  engine  serial  number.  These  times  are 
“not”  statistically  independent  because  of  the  manner  in  which  aircraft  engines  are  repaired. 
More  often  than  not,  they  are  not  totally  overhauled  before  being  reinstalled.  Only  minor  repair 
is  made  on  some  occasions.  The  operating  time  between  the  first  and  second  failure  of  an 
engine  is  dependent  on  the  operating  time  to  first  failure  and  the  extent  of  the  repair.  Engines 
with  long  first  failure  times  and  minor  first  repair  will  contain  many  used  components  with  large 
accumulated  operating  hours  when  reinstalled.  They  will  likely  have  shorter  operating  times  to 
the  next  failure  than  engines  with  small  operating  times  to  first  failure.  Consequently,  the 
operating  times  between  engine  failures,  for  the  same  engine,  are  not  statistically  independent. 
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B.  FrmNG  THE  PROBABH^ITY  DISTRIBUTION  OF  OPERATING  TIME 

BETWEEN  FAILURES  FOR  TF-34  ENGINES 

1.  Weibull  Distribution 

Aircraft  engines  wear  out.  Consequently,  any  probability  distribution  used  to  model 
their  time  to  failure  should  reflect  wear-out  characteristics.  The  Weibull  distribution  has  these 
properties.  It  was  used  by  Caudill  [Ref  1],  Malsbury  [Ref  2],  and  Mlikan  [Ref  3]  to  model 
the  distribution  of  operating  times  between  failures  for  TF-34  engines.  They  used  a  modified 
method  of  moments  procedure  to  compute  the  parameters  of  the  fitted  Weibull  distribution. 
This  method  used  special  properties  of  the  Weibull  distribution  and  the  exponential  distribution; 
the  sample  mean  and  the  sample  standard  deviation  were  used  to  estimate  the  parameters  of  the 
fitted  Weibull  distribution.  It  was  chosen  because  it  was  known  to  be  less  sensitive  to  large 
positive  or  negative  correlations  (dependency)  in  the  data  than  the  more  standard  methods  of 
maximum  likelihood  estimation  when  statistical  independence  is  assumed. 

2.  Research  Question 

The  analysis  performed  in  this  thesis  provides  a  comparison  of  the  accuracies  of  the 
maximum  likelihood  estimation  procedure  (assuming  independence)  and  the  modified  method 
of  moments  procedure  used  by  Caudill  for  fitting  a  Weibull  distribution  to  the  probability 
distribution  of  operating  times  between  engine  failures.  These  comparisons  are  made  for  times 
that  are  randomly  generated  in  a  manner  that  resembles  the  way  in  which  engine  failure  times 
occur  in  the  real  world.  These  randomly  generated  operating  times  between  failures  for  the 
same  engine  will  not  be  statistically  independent.  Once  this  comparison  was  completed,  it  was 
extended  to  real  data  using  TF-34  engine  operating  times  between  failures  from  the  NALDA 
data  base.  [Ref  1] 
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n.  GENERATION  OF  FAILURE  TIMES 


A.  PROCEDURE 

To  generate  the  failure  times,  a  program  in  S-Plus  (See  Appendix  A)  was  developed 
with  the  following  properties. 

1.  Each  engine  has  3  components  (subsystems). 

2.  The  failure  time  distributions  of  each  of  these  components  are  Weibull  distributions 
with  fixed  parameters  determined  a-priori  and  based  on  experience. 

3.  The  data  generated  is  a  Type  n  Censored  Data  Set,  and  a  total  of  3  failure  times 
are  gathered  for  each  engine. 

With  these  criteria,  the  program  generates  the  failure  times  in  the  following  way: 

Step  1:  Three  Weibull  random  numbers  are  generated.  They  represent  the  first  failure 
time  that  each  one  of  the  three  components  could  have.  These  values  constitute  the 
vector  faultsl.  Three  additional  Weibull  random  numbers  are  generated.  They 
represent  the  second  failure  time  that  each  one  of  the  three  components  could  have. 
These  values  constitute  the  vector  faultsHDisQt  additional  Weibull  random  numbers 
are  generated.  They  represent  the  third  failure  time  that  each  one  of  the  three 
components  could  have.  These  values  constitute  the  vector  faultsS. 

Step  2:  The  minimum  value  of  the  vector  faultsl  is  selected,  this  value  is  the  time  of 
the  first  failure  of  the  engine.  The  program  creates  a  new  vector  newfaill  that  is  the 
difference  between  the  vector  faultsl  and  the  minimum  value  selected;  it  is  clear  that 
this  new  vector  has  a  component  that  is  0.  This  component  is  replaced  in  the  vector 
newfaill  by  the  element  of  the  vector  faultsl  that  has  the  same  index  as  the  first  failure. 

Step  3:  At  this  point  the  engine  has  two  components  that  have  accumulated  working 
hours  and  one  that  has  been  just  repaired.  The  vector  newfaill  has  two  elements  that 
have  the  remainder  of  the  theoretical  first  failure  for  these  components  and  one  element 
that  is  the  theoretical  second  failure  for  the  other  component.  So,  again,  the  minimum 
element  of  this  vector  is  selected  and  it  is  the  second  failure  time  of  the  engine. 
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Step  4:  The  program  creates  a  new  vector  called  newfailS  that  is  the  difference 
between  the  vector  newfaiU  and  this  minimum  value  selected;  it  is  clear  that  this  vector 
has  a  component  that  is  0. 

If  the  second  engine  failure  time  has  taken  place  in  the  same  component  that  had  the  first 
failure,  the  0  element  is  replaced  by  the  element  of  the  vector  faults3  whose  index  matches  that 
of  the  component  that  has  had  the  two  first  failures  or  by  the  element  of  the  vector  faults! 
whose  index  matches  with  the  component  that  had  the  second  failure.  Again,  there  is  an 
element  in  the  vector  faultsS  that  is  a  theoretical  failure  time  for  one  component  just  repaired, 
and  two  elements  that  have  the  remainder  of  theoretical  failures.  So,  once  again  the  minimum 
element  of  this  vector  will  be  the  value  of  the  third  engine  failure  time.  The  vector  enginefails 
gives  the  three  random  failure  times  of  the  engine. 

B.  GENERATION 

An  S-Plus  program  was  developed  to  generate,  in  the  bank  of  engines,  the  three  failure 
times  for  the  engine.  The  program  can  be  mn  for  different  numbers  of  engines.  In  this  thesis  the 
program  was  mn  for  100  engines  with  a  total  of  300  failure  times  and  for  20  engines  with  60 
failure  times.  This  was  done  to  determine  if  the  number  of  engines  has  an  apparent  effect  in  the 
results  of  the  study.  For  each  one  of  these  engine  sets,  30  replications  were  generated.  These 
replications  were  exported  fi'om  S-Plus  to  Excel  with  the  export  tool  via  the  S-Plus  file-Export 
menu  and  used  to  assess  averages  on  the  accuracy  of  the  two  methods  of  estimation  • 

The  first  replication  for  both  20  and  100  engines  can  be  seen  in  the  Appendk  A,  all  of 
the  other  replications  have  been  retained  on  diskettes  attached  to  the  master  copy  of  this  thesis. 
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ffl.  ESHMAHON  OF  PARAMETERS  WITH  THE  MAXIMUM  LIKELIHOOD 
AND  THE  MODIFIED  MOMENT  METHOD 


A.  PROCEDURE 

The  two  parameter  Weibull  distribution  analyzed  in  this  thesis  has  a  probability  density 
function  given  by: 

m  =  ,  (1) 

or  equivalently 

M)  =  p0-l‘tf-'e‘'^  ,  (2) 


where  6  =  :r  -  Customarily  3  is  called  the  shape  parameter  and  X  (or  0 )  is  called  the  scale 

A/ 

parameter. 

To  estimate  the  parameters,  the  capabilities  of  the  spread  sheet  program  Excel  were 
used.  The  failure  times  are  generated  in  S-Plus  and  exported  to  Excel,  where  they  are  displayed 
automatically  on  the  Sheet  1  Worksheet  of  Excel. 

The  Maximum  Likelihood  estimators  are  calculated  based  on  well  known  formulas 

[Ref  4].  Solve  numerically  for  3 ; 


Y,LnTi- 


n^Ti^LnTi 


/=1 


(3) 


Then: 
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1 


(4) 


(5) 


The  parameters  for  the  Modified  Method  of  Moments  were  calculated  based  on  the 
fact  that  if  r  has  a  Weibull  distribution  then  has  a  Exponential  distribution  with  a  scale 
parameter  of  XP .  The  following  demonstration  is  fi'om  [Ref  5]: 

Considering  the  Weibull  distribution  given  by: 

m = .  (6) 


Let  x  =  so  t-  , 


dt  1  ^-1 


(7) 


Then 


f^{x)  =  I  J|/(?(jc))  = 


P 


1  i-i  ^ 

P 


(8) 


8 


Therefore: 


X  ~  exp(/l^) . 


(9) 


It  is  well  known  that  the  mean  and  standard  deviation  of  the  exponential  distribution 
are  equal.  Therefore  their  statistical  estimates  i?  and  s  should  be  nearly  equal.  That  is 


n-ltt 


2  1  ^ 


n 


(10) 


1=1 


where  the  T.  are  the  times  between  engine  failures. 


The  solution  for  P  in  equation  (10)  will  be  denoted  by  P . 

f  1  " 

Since  E(Tf )  =  6  ~  ,  the  sample  mean,  — T  Tf  estimates  9  ^ ,  therefore: 

VA./ 


e= 


1  ” 


L«  r=l 


(11) 


This  is  the  same  estimate  obtained  for  the  MLE  of  0. 


B.  ESTIMATION 

To  carry  out  the  estimation  of  parameters  with  the  two  methods,  two  Macros  in  Visual 
Basic  were  developed.  One  Macro  for  20  engines  (See  Appendix  B),  that  uses  60  failure  times, 
and  the  other  for  100  engines  (See  Appendbc  C),  that  uses  300  failure  times. 

Basically,  these  two  Macros  are  similar  and  therefore  they,  and  procedures  for  using 
them,  can  be  similarly  described.  Once  the  data  is  exported  to  Excel,  it  appears  on  Sheetl.  At 
this  point  it  is  necessaiy  to  open  the  file  with  the  Macro,  called  Macro20.xls  for  the  sample  of 
20  engines  and  Macrol00.xls  for  the  sample  of  100  engines  respectively.  Next  the  window 
with  the  Macro  is  hidden  fi-om  the  same  Sheetl. Th&n  the  Macro  can  be  run.  To  run  the  Macro, 
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choose  Macro  from  the  Tools  menu.  From  Macro/Name  Reference  choose  the  name  of  the 
Macro  and  click  on  run.  The  Macro  does  all  the  previous  calculations  and  operations  needed 
to  estimate  the  parameters  and  to  entitle  the  entire  Worksheet.  After  these  calculations  the 
Macro  stops.  At  this  point  it  is  necessary  to  introduce  the  following  data: 

In  Cell  B6  introduce  the  value  previously  calculated  by  the  macro  for  the  initial 
value  of  the  shape  parameter  which  is  in  Cell  B4. 

This  initial  parameter  was  introduced  to  assure  the  convergence  with  the  Tool  Goal 
Seeker,  because  without  it,  a  solution  may  not  be  reached.  The  initial  value  is  given  by  [Ref 
6]: 


(■ 

6 

n  n 

'ZilnTif-iZlnTif  In 

>.  /=1  /=1 

n-l 

1 

2 


(12) 


Using  these  choices  for  the  initial  shape  parameter,  an  average  of  only  3.5  Newton 
iterations  were  needed  to  achieve  four-place  accuracy  in  the  solution  for  p . 

After  introducing  the  initial  value  given  in  Cell  B4  into  Cell  B6,  cell  C6  must  be 
highlighted.  Then  the  Tools  Menu  is  selected.  Goal  Seek  is  selected,  a  0  is  placed  in  the  “To 
Value,”  and  B6  is  entered  as  the  “Changing  Cell.” 

Equation  (3)  is  in  Cell  C6  and  the  Goal  Seeker  searches  for  the  value  of  P  which 
satisfies  that  equation.  There  were  no  non-convergence  problems  in  any  of  the  70 
replications  made  with  that  formula.  Once  the  solution  is  achieved,  the  final  value  for  the 
shape  parameter  appears  in  Cell  B6.  The  solutions  for  X  and  0  appear  in  Cell  B9 
(Equation  (4))  and  B1 1  (Equation  (5))  respectively.  These  solutions  for  X  and  0  are 
the  maximum  likelihood  estimates  for  these  parameters  of  the  Weibull  distribution  that  is 
used  to  model  the  probability  distribution  of  the  family  of  times  between  failures  on  all 
engines.  These  values  are  printed  on  the  left  side  of  Sheet2. 

On  the  right  side  of  Sheet  2  are  the  calculations  carried  out  by  the  Modified 
Moment  Method. 
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There  were  convergence  problems  in  calculating  P  using  Equation  (10).  Amazingly 
the  same  formula  may  converge  or  not  depending  on  the  terms  on  each  side.  Through  trial  and 
error  it  was  discovered  that  by  writing  the  formula  in  the  following  way  it  was  possible  to  get 
convergence  in  all  replications: 


It-,' 


(13) 


In  Cell  Q6  it  is  necessary  to  introduce  the  initial  shape  parameter  given  in  Cell  Q4 
(Equation  (12)).  Once  done,  highlight  Cell  P7  (Equation  (13)),  select  the  Tools  Menu,  select 
Goal  Seek  and  depending  on  the  number  of  engines,  enter  in  “To  Value”  7.811334  (20 
engines,  60  failure  times,  60/ )  or  17.349447  (100  engines,  300  failure  times,  300/-\/299  ). 
In  the  “Changing  Cell,  Q6,”  the  Goal  Seeker  looks  for  the  P  value  that  will  satisfy  Equation 
(13).  Cells  Q9  and  Qll  contain  the  values  for  the  scale  parameter  X  and  0  determined  by 
Equation  (1 1)  and  (5). 

A  total  of  30  replications  for  each  set  of  engines  were  done.  Afterwards,  for  the 
reasons  explained  later,  10  more  replications  were  done  using  real  data. 

Printouts  of  the  failure  times  and  parameter  calculations  for  the  20  and  100  engines  first 
replication  are  provided  in  Appendix  E. 

A  summary  of  the  parameters  obtained  is  shown  in  the  following  tables; 
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Table  1.  Twenty  Engines  =  60  Failure  Times 


20  Engines  = 

Maximum  Likelihood  Estimator 
Shape  Scale 


60  Failure  Times 

Modified  Moment  Method 
Shape  Scale 


11 

0.973000000 

251.705300000 

m 

1.105104390 

279.481900000 

E 

1.052707000 

242.343600000 

ri 

1.005354000 

282.577600000 

m 

1.211857000 

272.413900000 

1.523761000 

277.443800000 

ii 

1.053815000 

243.377800000 

1.015065000 

330.527000000 

m 

0.961137000 

251.822700000 

E 

1.108272000 

263.760200000 

H 

1.258861000 

317.378200000 

H 

1.358459000 

245.953800000 

E 

1.323305000 

281 .038800000 

E 

1.054978965 

223.276800000 

H 

1.228101693 

264.364500000 

1.161814103 

243.734600000 

1.355617000 

268.067900000 

1.308393000 

309.037200000 

E 

1.391697000 

284.440700000 

E 

1.141031467 

242.289400000 

1.336786000 

271.132800000 

E 

1.389460000 

244.898600000 

E 

1.200747000 

288.213200000 

E 

^  1.297427000 

293.585500000 

■  1.400818000 

310.907900000 

E 

1.231510000 

277.128900000 

E 

1.048038281 

277.331100000 

m 

1.093433000 

273.733700000 

E 

1.295133208 

246.789500000 

E 

1.224688831 

303.055100000 

Average  1.203679065  272.060400000 
Std.  Dev.  0. 1 4588671 6  25.376732459 


0.985725000 

253.139700000 

1.129925000 

281 .796500000 

1.036553000 

240.792500000 

1.079900000 

291 .030800000 

1.102670000 

263.279100000 

1.442500000 

273.090800000 

1.309674000 

262.946000000 

1.006751000 

329.425200000 

1.078819000 

264.611000000 

1.084296000 

261.569300000 

1.395249000 

328.224900000 

1.211490000 

237.176300000 

1.461631834 

289.757400000 

1.249979811 

236.965600000 

1.182801000 

260.845800000 

1.239576000 

249.537300000 

1.367488000 

268.798900000 

1.357252000 

312.609300000 

1.406588000 

285.359900000 

1.324936416 

255.330700000 

1.223524000 

263.236700000 

1.280865198 

238.618200000 

1.263844000 

'  293.349300000 

1.260780021 

290.891900000 

1.325062000 

mmmm\ 

1.320218000 

283.439100000 

1.043436000 

276.830600000 

1.248177000 

286.965800000 

1.240951783 

243.300600000 

1.402646081 

317.067700000 

1 .235443638  274.848080000 
0.1 33749260  25.937151 577 
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Replication  Number 


Table  2.  One-Hundred  Engines  =  300  Failure  Times 


100  Engines  =  300  Failure  Times 


Maximum  Likelihood  Estimator 


Modified  Moment  Method 


1.261028000 

292.556900000 ! 

M 

1.192249000 

267.935700000 

m 

1.191393000 

301.254500000 

m 

1.100167000 

265.867700000 

m 

1.132621000 

273.303500000 

B 

1.220806000 

282.032200000 

■ 

1.206135000 

267.017000000 

1 .254734000 

291.692500000 

M 

1 .284450000 

285.463200000 

1 .277785000 

285.850800000 

1 .274025000 

263.121900000 

1.194163000 


1.216995000 


IH335!SB!»lili 


1.416971000 


1 .264662000 


1.342221000 


1 .229806000 


1.199120000 


1.189326000 


1.257310000 


267.935700000 


271.996400000 


272.494100000 


293.370600000 


255.640200000 


294.799600000 


268.682800000 


279.712400000 


287.109500000 


267.635600000 


288.996500000 


263.550100000 


261 .778100000 


256.670600000 


1 .127942000  279.203900000 


Average  1 .226599300  275.1 85723333 
Std.  Dev.  0.063780678  12.379509005 


1.288232000 

294.659500000 1 

1.182793000 

1.250620000 

272.423600000 

1.243339000 

305.757700000 

1.187671000 

272.564100000 

1.205131000 

279.422100000 

1.334481000 

290.504100000 

1 .246956000 

270.084700000 

1.271828000 

1.292318000 

286.039200000 

1 .387920000 

293.508300000 

1.341512000 

267.591700000 

1 .250620000 

272.423600000 

1 .228539000 

274.677200000 

1.231497000 

273.612100000 

1 .200586000 

296.122500000 

1 .259945000 

254.068200000 

1 .422087000 

295.112700000 

1 .379585000 

276.306900000 

1 .448534000 

286.372700000 

1 .269967000 

290.254800000 

1.281756000 

273.804000000 

1.207016000 

290.494900000 

1.346121000 

273.747200000 

1.201983000 

256.990000000 

1 .249969000 

256.167900000 

1.311836000 

273.423600000 

1 .237968000 

278.973700000 

1.219361 000 

263.912300000 

1.173691 000 

283.295300000 

1.271795400  278.553426667 
0.07051 61 79  1 3.009908323 
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IV.  ACCURACY  ANALYSIS 


A.  PROCEDURE 

The  accuracy  of  each  method  was  measured  using  the  Kolmogorov-Smimov  goodness 
of  fit  Test  and  the  S-Plus  software  package.  The  null  hypothesis  in  S-Plus  is  “the  True  cdf 
(empirical)  equals  the  Weibull  distribution”.  This  test  has  the  advantage  that  it  uses  continuous 
populations.  Therefore  it  is  not  necessary  to  group  the  data. 

This  test  has  the  disadvantage  that  when  distribution  parameters  are  estimated  fi'om  the 
same  sample  of  data,  the  test  produces  a  conservative  result;  that  is,  the  probability  of  a  Type  I 
error  will  be  smaller  than  specified  [Ref  6].  This  drawback  has  only  a  relative  importance  in  the 
present  study.  The  K-S  Test  p-value  is  a  conservative  result  whether  the  parameters  are 
estimated  by  the  Maximum  Likelihood  Method  or  the  Modified  Moment  Method. 

A  S-Plus  firnction  called  Tocompare  (See  Appendbc  E)  was  created  to  calculate  the  K- 
S  statistics  for  the  set  of  parameters  calculated  in  each  replication  using  both  methods.  This  S- 
Plus  function  computes  the  K-S  Test  Statistic  (Supremum)  and  the  p-value  for  each  method, 
and  it  also  plots  the  empirical  and  hypothesized  distributions.  In  Appendbc  E  the  reader  will  find 
the  output  of  the  function  for  the  first  replication  of  each  set  of  parameters. 

B.  ACCURACY 

The  calculations  for  all  replications  and  the  values  obtained  for  each  set  of  engines  are 
presented  in  Tables  3  and  4.  The  letter  “  B”  marks  the  method  that  achieved  a  better  fit,  using 
as  reference  the  K-S  statistic;  that  is,  the  method  with  the  smaller  K-S  statistic  is  the  more 
accurate  procedure  for  that  replication. 

For  the  set  of  20  engines,  the  MLE  method  gave  a  better  fit  in  17  of  the  30  replications. 
For  the  set  of  100  engines,  the  MLE  method  gave  a  better  fit  in  20  of  the  30  replications. 
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Table  3.  Better  Fit  Twenty  Engines 


20  Engines  =  60  Failures  times 

M.L.E.  M.M.M. 

shape  scale  k-s  p-value  shape  scale  k-s  p-value 


1  j  0.973000000 
1.105104390 
1.052707000 
1.005354000 
1.211857000 
1.523761000 
1.053815000 
1.015065000 
0.961137000 
1.108272000 
1.258861000 
1.358459000 
1.323305000 
1.054978965 
1.228101693 
1.161814103 
1.355617000 
1.308393000 
1.391697000 
1.141031467 
1.336786000 
1.389460000 
1.200747000 
24]  1.297427000 
1.400818000 
1.231510000 
2?j  1.048038281 
1.093433000 
1.295133208 
1.224688831 


251.705300000 

279.481900000 

242.343600000 

282.577600000 

272.413900000 

277.443800000 

243.377800000 

330.527000000 

251.822700000 

263.760200000 

317.378200000 

245.953800000 

281.038800000 

223.276800000 

264.364500000 

243.734600000 

268.067900000 

309.037200000 

284.440700000 

242.289400000 

271.132800000 

244.898600000 

288.213200000 

293.585500000 

310.907900000 

277.128900000 

277.331100000 

273.733700000 

246.789500000 

303.055100000 


0.054600000 

0.088300000 

0.076300000 

0.061700000 

0.101800000 

0.086700000 

0.091500000 

0.052700000 

0.088000000 

0.064500000 

0.101500000 

0.104300000 

0.075700000 

0.088000000 

0.070300000 

0.078800000 

0.070200000 

0.049400000 

0.083700000 

0.091800000 

0.107000000 

0.085000000 

0.083700000 

0.058400000 

0.071100000 

0.062500000 

0.078500000 

0.078000000 

0.055600000 

0.083900000 


0.994000000 

0.737300000 

0.876000000 

0.976300000 

0.562800000 

0.757300000 

0.696500000 

0.996200000 

0.742000000 

0.964400000 

0.566300000 

0.531000000 

0.881800000 

0.741900000 

0.928000000 

0.850400000 

0.929200000 

0.998600000 

0.794200000 

0.692600000 

0.498400000 

0.778500000 

0.795100000 

0.986700000 

0.921800000 

0.973200000 

0.853300000 

0.859200000 

0.992500000 

0.792200000 


mean 

mean 

mean 

mean 

1.203679065 

272.060400000 

0.078116667 

0.822256667  17 

stdev 

stdev 

stdev 

stdev 

0.145886716 

25.376732459 

0.015610639 

0.146004782 

1 0.985725000 
1.129925000 
1.036553000 
bI  1.079900000 
*  1.102670000 
1.442500000 
B  1.309674000 
B  1.006751000 
B  1.078819000 
1.084296000 
B  1.395249000 
B  1.211490000 
B  1.461631834 
1.249979811 
1.182801000 
B  1.239576000 
B  1.367488000 
B  1.357252000 
1.406588000 
Bl  1.324936416 
1.223524000 
B  1.280865198 
B  1.263844000 
B  1.260780021 
1.325062000 
1.320218000 
1.043436000 
Bl  1.248177000 
Bj  1.240951783 
bI  1.402646081 


mean 

1.235443638 

stdev 

0.133749260 


253.139700000 

281.796500000 

240.792500000 

291.030800000 

263.279100000 

273.090800000 

262.946000000 

329.425200000 

264.611000000 

261.569300000 

328.224900000 

237.176300000 

289.757400000 

236.965600000 

260.845800000 

249.537300000 

268.798900000 

312.609300000 

285.359900000 

255.330700000 

263.236700000 

238.618200000 

293.349300000 

290.891900000 

305.455500000 

283.439100000 

276.830600000 

286.965800000 

243.300600000 

317.067700000 


mean 

274.848080000 

stdev 

25.937151577 


0.050900000 

0.081000000 

0.073000000 

0.086100000 

0.097300000 

0.081000000 

0.115700000 

0.053100000 

0.098600000 

0.061500000 

0.126600000 

0.121500000 

0.110800000 

0.087200000 

0.065700000 

0.101500000 

0.072500000 

0.054400000 

0.083400000 

0.116800000 

0.077500000 

0.092700000 

0.096600000 

0.059800000 

0.071100000 

0.057400000 

0.077100000 

0.084100000 

0.067600000 

0.095900000 


mean 

0.083946667 

stdev 

0.020788358 


0.997700000  B| 

0.826600000  B 

0.906600000  B 

0.765600000 

0.621600000  B 

0.826400000  B 

0.398400000 

0.995900000 

0.604700000 

0.977100000  B 

0.290900000 

0.338000000 

0.452600000 

0.751700000  B 

0.958200000  B 

0.567200000 

0.910600000 

0.994300000 

0.798500000  B 

0.386000000 

0.863800000  B 

0.681600000 

0.629800000 

0.982700000 

0.922100000  B 

0.989000000  B 

0.867800000  B 

0.790300000 

0.947000000 

0.639300000 


mean 

0.756066667  13 
stdev 
0.214131410 
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Table  4.  Better  Fit  One-Hundred  Engines 


100  Engines  =  300  Failures  times 

M.L.E.  M.M.M. 

k-s  p-value  shape  scale  k-s  p>va!ue 


111.261028000 

1.147177000 

1.192249000 

1.191393000 

1.100167000 

1.132621000 

1.220806000 

1,206135000 

1.254734000 

1.284450000 

1.277785000 

1.274025000 

1.192249000 

1.194163000 

1.216995000 

1.169718000 

1.283386000 

1.416971000 

1.264662000 

1.342221000 

1.229806000 

1.199120000 

1.189326000 

1.199541000 

1.268509000 

1.257310000 

1.274659000 

1.208075000 

1.220756000 

1.127942000 


292.556900000 

258.382300000 

267.935700000 

301.254500000 

265.867700000 

273.303500000 

282.032200000 

267.017000000 

291.692500000 

285.463200000 

285.850800000 

263.121900000 

267.935700000 

271.996400000 

272.494100000 

293.370600000 

255.640200000 

294.799600000 

268.682800000 

279.712400000 

287.109500000 

267.635600000 

288.996500000 

263.550100000 

261.778100000 

256.670600000 

270.850700000 

276.652800000 

264.013900000 

279.203900000 


0.031900000 

0.027400000 

0.047000000 

0.055800000 

0.060600000 

0.043200000 

0.034300000 

0.037900000 

0.021600000 

0.022700000 

0.044600000 

0.047400000 

0.029400000 

0.039600000 

0.031400000 

0.034800000 

0.031800000 

0.039800000 

0.053800000 

0.039700000 

0.041300000 

0.036000000 

0.032800000 

0.032000000 

0.028200000 

0.026400000 

0.038700000 

0.035200000 

0.033800000 

0.035700000 


0.920600000  B| 
0.977600000  B 
0.522300000  B 
0.306900000 
0.220000000 
0.630900000 
0.873100000  B 
0.782600000 
0.999000000  B 
0.997900000 
0.589000000  B 
0.511700000  B 
0.957700000  B 
0.735300000  B 
0.928900000 
0.861700000  B 
0.922000000  B 
0.729300000 
0.349300000  B 
0.730800000  B 
0.685900000 
0.830800000  B 
0.903200000  B 
0.919200000  B 
0.971000000  B 
0.984800000  B 
0.759500000  B 
0.851300000 
0.882800000  B 
0.838800000 


1.288232000 

1.182793000 

1.250620000 

1.243339000 

1.187671000 

1.205131000 

1.334481000 

1.246956000 

1.271828000 

1.292318000 

1.387920000 

1.341512000 

1.250620000 

1.228539000 

1.231497000 

1.200586000 

1.259945000 

1.422087000 

1.379585000 

1.448534000 

1.269967000 

1.281756000 

1.207016000 

1.346121000 

1.201983000 

1.249969000 

1.311836000 

1.237968000 

1.219361000 

1.173691000 


294.659500000 

261.246400000 

272.423600000 

305.757700000 

272.564100000 

279.422100000 

290.504100000 

270.084700000 

293.041800000 

286.039200000 

293.508300000 

267.591700000 

272.423600000 
274.677200000 
273.612100000 
296.122500000 
254.068200000 
295.112700000 
276.306900000 
286.372700000 
290.254800000 
273.804000000 
290.494900000 
273.747200000 
256.990000000 
256.167900000 

273.423600000 
278.973700000 
263.912300000 
283.295300000 


0.038200000 

0.032700000 

0.051000000 

0.046600000 

0.045600000 

0.042500000 

0.048900000 

0.036300000 

0.025500000 

0.020800000 

0.058300000 

0.048900000 

0.035300000 

0.048600000 

0.029400000 

0.041300000 

0.033000000 

0.038700000 

0.065700000 

0.067500000 

0.040300000 

0.044400000 

0.033600000 

0.048900000 

0.041500000 

0.027700000 

0.048400000 

0.028200000 

0.034200000 

0.031400000 


0.774000000  I 

0.905100000 

0.416900000 

0.532800000  B 

0.561900000  B 

0.651400000  B 

0.469300000 

0.823000000  B 

0.989900000 

0.999500000  B 

0.260200000 

0.469900000 

0.848500000 

0.477100000 

0.957500000  B 

0.685300000 

0.898900000 

0.759500000  B 

0.150500000 

0.129900000 

0.713200000  B 

0.596300000 

0.886800000 

0.469200000 

0.680900000 

0.975800000 

0.482200000 

0.971300000  B 

0.874100000 

0.928300000  B 


mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

1.226599300 

275.185723333 

0.037160000 

0.772463333  20 

1.271795400 

278.553426667 

0.041113333 

0.677973333  10 

standev 

standev 

standev 

standev 

standev 

standev 

standev 

standev 

0.063780678 

12.379509005 

0.009086679 

0.209187169 

0.070516179 

13.009908323 

0.010947594 

0.247541499 
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If  the  MLE  and  MMM  methods  are  equally  good  then  one  would  expect  that  the 
probability  the  MLE  gives  a  better  fit  (  using  the  K-S  statistic)  should  be  0.5.  Let  tt  represent 
the  probability  the  MLE  gives  a  better  fit;  the  above  sets  of  data  can  then  be  used  to  test  Ho :  n 
=  0.5  versus  Ha ;  7t  0.5.  For  the  set  of  20  engines,  17  successes  were  observed  out  of  30 
replications  ( MLE  gave  a  smaller  K-S  value  than  MMM);  the  p-value  for  testing  Ho :  tt  =  0.5 
is  0.5847  so  Ho  would  be  accepted  for  any  a  <  0.5847.  For  the  set  of  100  engines,  20  of  the 
30  replications  gave  successes;  in  this  case  the  p-value  is  0.0987  and  Ho  is  rejected  for  any  a 
>  0.0987.  Certainly  this  comparison  of  MLE  and  MMM,  using  the  binomial  test,  shows  no  real 
difference  for  20  engines  and  marginally  indicates  superiority  of  MLE  for  100  engines. 

If  the  K-S  statistic  mean  and  standard  deviation  are  analyzed  it  is  clear  that  the  MLE. 
gives  the  smaller  mean  and  standard  deviation  of  the  two  cases. 

For  this  data,  it  seems  that  despite  the  fact  that  there  is  some  dependence  between  the 
failure  times,  the  MMM  does  not  give  a  better  fit  than  the  MLE.  This  is  true  at  least  for  the 
failure  times  generated  Avith  the  simulation  program  used  in  this  study.  This  is  likely  due  to  the 
fact  that  the  failure  times  for  different  engines  (among  engines)  are  independent. 

At  this  point  real  failure  data  was  examined.  This  real  data  is  from  [Ref  1],  and  from 
the  cleansed  data  collected  in  Appendix  B  (Five-Year  Time  Between  Failures  Data  Base).  This 
data  was  cleansed,  removing  zeros  as  failure  times.  Ten  samples  of  data  were  taken  fi'om  the 
data  set  of  operating  times  between  failures  on  TF  34  aircraft  engines.  Five  samples  were  of  60 
failures  times  and  five  samples  were  of 300  failure  times. 

The  same  procedure  to  estimate  the  parameters  and  analyze  the  accuracy  of  the  results 
were  applied.  The  results  are  shown  in  Tables  5  and  6. 

With  the  real  data  the  result  is  the  same,  a  better  fit  is  obtained  more  often  with  the 
MLE  method  than  with  the  MMM  method.  In  the  60  failure  times  case  the  MLE  method  was 
the  better  procedure  in  four  of  the  five  samples,  and  in  the  300  failure  times  case  it  was  better 
three  of  the  five  samples  using  the  K-S  statistic  as  measure  of  accuracy. 
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It  is  important  to  point  out  that  in  these  real  data,  the  engine  structure  is  unknown. 
There  is  no  information  about  how  many  elements  the  engine  has,  or  of  course  the  way  that 
these  components  failed. 

Table  5.  Estimate  Parameters  and  Analyze  Accuracy  Twenty  Engines 


A  w 

.S  a> 

s  s 
■■3  s 
g-z 


20  Engines  =  60  Failures  times 

M.LE.  M.M.M. 

shape  scale  k-s  p-value  shape  scale  k-s  p-value 


1.554927  651.243400  0.103000  0.547400  B 

1.435165  710.748800  0.079700  0.840000  B 

1.531311  603.432400  0.072600  0.909800 

1.436419  609.475100  0.103800  0.537100  B 

1.535733  579.331200  0.135400  0.221300  B 

1.732451  669.495500  0.103200  0.544600 
1.735014  748.435600  0.130200  0.261300 
1.635563  613.677100  0.051800  0.997100  B 
1.485100  615.440400  0.114900  0.406200 
1.490634  574.988000  0.145900  0.155400 

mean  mean  mean  mean 

1.414134  642.431620  0.048360  0.511380  4 
stdev  stdev  stdev  stdev 

0.045711  14.890226  0.008406  0.185231 

mean  mean  mean  mean 

1.587754  663.677300  0.057660  0.289300  1 
stdev  stdev  stdev  stdev 

0.023623  16.834529  0.006280  0.138620 

C  ANALYSIS  OF  THE  RESULTS 

The  author  expected  “  a-priori”  a  better  fit  with  the  MMM  than  with  the  MLE,  based 
on  the  fact  that  the  generated  operating  times  between  failures  for  the  same  engine  are  not 
statistically  independent.  But  the  resulting  fit  for  the  data  was  better  with  the  MLE  more  often 
than  with  the  MMM.  This  result  is  further  analyzed  here,  by  examining  how  the  dependence 
between  failure  times  affect  each  method. 

To  assess  dependence  between  failure  times,  the  correlation  and  the  grade  of 
dependence  between  the  failure  times  was  calculated  for  each  replication.  The  correlation 
coefficients  were  calculated  and  added  at  the  end  of  each  replication  (See  Appendk  A).  In 
addition  the  Kendall-Tau  independence  test  was  performed.  The  null  hypothesis  in  S-Plus  for 
this  test  is  “  Tau  is  equal  to  0”  that  is  the  two  samples  are  independent.  This  analysis  was 
performed  for  all  possible  pairwise  combinations  of  the  failure  times  of  an  engine;  namely  first- 
second,  first-third,  second-third  correlations,  tau  coefficients  and  p-values  were  calculated.  For 
the  Kendall-Tau  test  the  minimum,  maximum  and  mean  p-values  of  those  pairwise 
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combinations  of  each  replication  were  added  to  the  tables.  The  results  are  shown  in  Tables  7 
through  10. 

Table  6.  Estimate  Parameters  and  Analyze  Accuracy  One-Hundred  Engines 


100  Engines  =  300  Failures  times 

M.LE.  M.M.M. 

Shape  scale  k-s  p-value  shape  scale  k-s  p-value 

©  5;  1  1.462681  618.916600  0.047200  0.515100  1.579829  632.033000  0.045700  0.559000  B 

'S’S  2  1.439461  663.159800  0.045200  0.571700  B  1.586168  681.405000  0.063000  0.184300 

.a  I  3  1.366409  641.062700  0.040500  0.709200  B  1.564422  666.063600  0.057300  0.279000 

g-Z  4  1.351920  636.911100  0.064600  0.163100  1.575530  665.621100  0.061900  0.200100  B 

0^  5  1.450201  652.107900  0.044300  0.597800  B  1.632821  673.263800  0.060400  0.224100 

mean  mean  mean  mean  mean  mean  mean  mean 

1.414134  642.431620  0.048360  0.511380  3  1.587754  663.677300  0.057660  0.289300  2 
stdev  stdev  stdev  stdev  stdev  stdev  stdev  stdev 

0.045711  14.890226  0.008406  0.185231  0.023623  16.834529  0.006280  0.138620 


From  those  tables  may  be  deduced  that  with  the  set  of  20  engines  the  correlation  is 
clearly  larger  than  with  100  engines.  When  the  independence  is  tested  with  the  Kendall-Tau 
test,  the  values  of  the  tau’s  coeflScients  are  similar  to  the  correlation  coefiBcients,  and  have 
larger  values  for  the  set  of  20  engines  than  the  set  of  100  engines.  In  general  for  this  set  of  data, 
the  correlation  coefficients  are  small  and  the  p-values  in  the  Kendall-Tau  are  not  significant 
enough  to  reject  the  null  hypothesis  of  independence  at  the  usual  values.  That  means  the 
dependence  obtained  in  the  simulation  is  small,  specially  with  the  set  of  100  engines;  this  may 
be  unfavorable  to  the  performance  of  the  MMM. 

To  summarize  the  results.  Figures  1  through  4  present  the  min  and  max  Kendall-Tau 
test  p-value  for  each  replication  for  the  method  that  gives  the  best  fit. 

From  these  figures  and  specially  fi-om  the  max  p-values,  it  may  be  deduced  that,  for  this 
data,  in  general  the  larger  the  max  p-value  (more  independence)  the  better  the  fit  with  the  MLE 
method.  Conversely  the  smaller  the  max  p-value  (less  independence)  the  better  the  fit  with  the 
MMM  method. 


Table  7.  Correlation  -  20  Engines 


Correlation  Coefficients  for  20  enaines 


rep16  Failurel  Fi 

Failurel  1  _ 

Failure2  0.381725000  _ 

Failures  O.i 


repi 

Failurel 

Failures 

Failures 

r 

Failurel 

1 

Failure2 

-0.099201099 

1 

Failures 

-0.325682845 

0.177625785 

1 

Failurel 

Failurel 
Failure2  >0,292364468 
Failures  -0.171746200 


Failurel 

1 

-0.125786346 

-0.240386929 


Failurel 

1 

-0.049045517 

0.141463053 


-0.146477812 


Failure2 


rep7 

Failurel 

Failures 

Failures 

■ 

Failurel 

1 

1 

Failures 

-0,334293594 

1 

1 

Failures 

1 

■ 

rep8 

Failurel  Failures 

Failures 

Failurel 

1 

Failures 

-0.362844249  1 

Failures 

-0. 15331 6474  -0.01 4631447 

1 

0.138519930 

-0.383549222 


Failurel 


0.096601004 


Failure2 


1 

-0.205353856 

1 

-0.202946234 

-0.062062633 

1 

Failurel  Failure2 

Failurel _ 1 _ 

Failure2  -0.431010000 _ 1 

Failures  -0.284860000  - 


Failurel  Fallure2 

Failurel _ 1 _ 

Failure2  -0.075920000 _ 1 

Failures!  0.111095000  -0.381490000 


Failure2  -0.019240000 
Failures  0.089556000  I  0.366561000 


Failurel 

Failurel _ 1 _ 

Failure2  -0.267500000 
Failures  0.140151000 


Falurel 

Failurel  1 

Failure2  0.368302000 
Failures  -0.205760000 


repSO 

Failurel 

Failures 

Failures 

Failurel 

1 

Failures 

-0.183030000 

1 

Failures 

-0,007070000 

-0.235020000 

1 

repSO 

Failurel 

Failures 

Failures 

Failurel 

1 

Failures 

0.011941000 

1 

Failures 

0.144591000 

0.137677000 

1 
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Table  8.  Correlation  -  100  Engines 


Correlation  Coefficients  for  100  engines 


repi 

Faiiirel 

Failire2 

Failures 

rep16 

Faiiirel 

FailireS 

FailireS 

Failurat 

1 

Failurel 

1 

Failure^ 

0.012345873 

1 

FailireS 

-0.343673022 

1 

Failures 

1 

FailireS 

-0.065777779 

-0.047403929 

1 

rap2 

Faiiirel 

Failire2 

FailireS 

rep17 

Failurel 

FailireS 

FailireS 

Failiral 

1 

Faiiirel 

1 

Failira2 

-0.053286726 

1 

FailireS 

-0.095841896 

1 

Failures 

-0.044274299 

1 

FailireS 

0.153445659 

-0.210881943 

1 

rap3 

Failurel 

Failire2 

FailireS 

repIS 

Failurel 

FailireS 

FailireS 

Failurel 

1 

Failurel 

1 

Failira2 

-0.250749723 

1 

FailireS 

0.052473163 

1 

FailiraS 

-0.094868699 

1 

FailireS 

0.041985512 

0.022318537 

1 

rap4 

Faiiirel 

Failures 

Failures 

rep19 

Failurel 

Failures 

Failures 

Failurel 

\mmam 

\mmm 

Faiiirel 

1 

Failure2 

1  -0.218109532 

1  1 

FaliuraS 

-0,062164561 

1 

Failures 

1 

Failures 

-0.1608692D1 

1 

rep5 

Failurel 

Failire2 

Failures 

repSO 

Failurel 

Failures 

FailireS 

Failtrel 

■■■1 

Failurel 

1 

Failure2 

1  -0.150008998 

1  1 

FailireS 

-0.083449627 

1 

FailireS 

1 

FailireS 

0.055670639 

1 

reps 

Faiiirel 

FailireS 

Failures 

repS1 

Faiiirel 

Failures 

FailireS 

Faiturel 

l■■HI 

Faiiirel 

1 

Failure2 

-0.079147073 

1 

Failures 

-0.110695458 

1 

FailiraS 

1 

FtulireS 

-0.170345154 

1 

rep7 

Faiiirel 

FailireS 

Failures 

repSS 

Failurel 

Failures 

Failures 

Failurel 

1 

Failurel 

1 

Failire2 

-0.065957419 

1 

Failures 

-0.180701469 

1 

FailiraS 

-0.065309916 

1 

FailireS 

0.156143924 

1 

rap8 

Faiiirel 

FailireS 

Failures 

repSS 

Faiiirel 

Failures 

Failures 

Faiiirel 

1 

Faiiirel 

1 

Failure2 

-0.087122634 

1 

FailireS 

-0.073156230 

1 

Failures 

0.009267376 

-0.198193315 

1 

Failures 

-0.171758743 

-0.083586000 

1 

rep9 

Failurel 

F»lure2 

FailireS 

rep24 

Failurel 

Failures 

Failures 

Faiiirel 

1 

Failurel 

1 

Fatlire2 

-0.024708152 

1 

FailireS 

-0.063590173 

1 

FailireS 

-0.088238686 

-0.153996216 

1 

Failures 

-0.106654719 

-0.162453112 

1 

rep  10 

Failurel 

FailireS 

Failures 

repSS 

Failurel 

Failures 

Failu-eS 

Failurel 

Failurel 

1 

Failire2 

-0.085599667 

1 

FailireS 

-0.084052612 

1 

FailireS 

-0.208242972 

-0.04475226 

1 

Failures 

0.049055891 

-0.063152872 

1 

repll 

Faiiirel 

FailireS 

FailireS 

repS6 

Failurel 

Failures 

Failures 

Failurel 

1 

Failurel 

1 

Failure2 

-0.030401682 

1 

FailireS 

-0.066404142 

1 

Failures 

0.019141802 

0.012598681 

1 

Failures 

0.099904638 

-0.219739207 

1 

rep12 

Faiiirel 

Failures 

FailireS 

rep27 

Faiiirel 

FailireS 

FailireS 

Faiiirel 

1 

Faiiirel 

1 

Failure2 

-0.243204236 

1 

FailireS 

-0.056026419 

1 

Failures 

0.016811427 

-0.06918138 

1 

FailireS 

-0.077399648 

-0.002855117 

1 

rep13 

Faiiirel 

FailireS 

FalureS 

repSS 

Failurel 

Failures 

Failu^eS 

Faiiirel 

1 

Failurel 

1 

Failire2 

-0.250749723 

1 

Failures 

-0.081782938 

1 

FailireS 

0.062100274 

-0.094868699 

1 

FailireS 

-0.061450025 

1 

rep  14 

Failurel 

FulireS 

FailireS 

repS9 

Failurel 

FailireS 

Failures 

Failurel 

BHOHII 

Failurel 

1 

Failure2 

-0.064039892 

1 

Failures 

-0.243170112 

1 

FailireS 

-0.077445912 

-0.077941408 

1 

Failures 

0.093441256 

-0.10701961 

1 

repIS 

Faiiirel 

FailireS 

FulireS 

repSO 

Fzulirel 

F»lureS 

Failures 

Faiiirel 

1 

Faiiirel 

1 

FeulureS 

-0.098178961 

1 

FailireS 

-0.23939123 

1 

FailireS 

EBSII 

1 

Failures 

-0.020603539 

-0.257213197 

1 
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Replication  Number 


Table  9.  Kendall  -Tau  Test  -  20  Engines 


tau12 

pvalue 

tau13 

pvalue 

tau23 

pvalue 

mean 

min 

max 

best  method 

ii 

-0.084 

0.604 

-0.242 

0.136 

0.126 

0.436 

0.392 

0.136 

0.604 

m.m.m 

ii 

-0.21 1 

0.194 

0.221 

0.173 

-0.062 

0.697 

0.355 

0.173 

0.697 

m.m.m 

M 

-0.116 

0.475 

-0.158 

0.330 

-0.095 

0.559 

0.455 

0.330 

0.559 

m.m.m 

ii 

-0.179 

0.270 

0.105 

0.516 

0.126 

0.436 

0.408 

0.270 

0.516 

m.I.e. 

ill 

-0.126 

0.436 

-0.295 

0.069 

-0.074 

0.650 

0.385 

0.069 

0.650 

m.m.m 

6 

0.105 

0.516 

0.105 

0.516 

0.221 

0.173 

0.402 

0.173 

0.516 

m.m.m 

7 

-0.232 

0.153 

-0.011 

0.948 

0.084 

0.604 

0.568 

0.153 

0.948 

m.I.e. 

8 

-0.221 

0.173 

-0.063 

0.697 

-0.168 

0.299 

0.390 

0.173 

0.697 

m.I.e. 

@ 

0.063 

0.697 

-0.295 

0.069 

0.179 

0.270 

0.345 

0.069 

0.697 

m.I.e. 

10 

-0,011 

0.948 

-0.274 

0.092 

-0.126 

0.436 

0.492 

0.092 

0.948 

m.m.m 

11 

-0.147 

0.364 

0.010 

0.948 

-0.063 

0.697 

0.670 

0.364 

0.948 

m.I.e. 

ii 

-0.116 

0.475 

-0.158 

0.330 

-0.095 

0.559 

0.455 

0.330 

0.559 

m.I.e. 

13 

-0.042 

0.795 

-0.179 

0.270 

0.147 

0.364 

0.476 

0.270 

0.795 

m.I.e. 

14 

-0.263 

0.105 

-0.179 

0.270 

-0.179 

0.270 

0.215 

0.105 

0.270 

m.m.m 

IS 

-0.168 

0.299 

0.232 

0.153 

-0.326 

0.044 

0.166 

0.044 

0.299 

m.m.m 

16 

0.189 

0.243 

0.221 

0.173 

-0.021 

0.897 

0.438 

0.173 

0.897 

m.I.e. 

17 

-0.126 

0.436 

-0.400 

0.014 

0.116 

0.475 

0.308 

0.014 

0.475 

m.I.e. 

16 

0.105 

0.516 

0.053 

0.746 

0.105 

0.746 

0.669 

0.516 

0.746 

m.I.e. 

n 

-0.126 

0.436 

-0.063 

0.697 

0.011 

0.948 

0.694 

0.436 

0.948 

m.m.m 

20 

-0.116 

0.475 

0.221 

0.173 

-0.200 

0.218 

0.289 

0.173 

0.475 

m.I.e. 

21 

-0.168 

0.218 

0.274 

0.092 

-0.095 

0.559 

0.289 

0.092 

0.559 

m.m.m 

22 

-0.21 1 

0.194 

-0.095 

0.559 

0.042 

0.795 

0.516 

0.194 

0.795 

m.I.e. 

26 

-0.189 

0.243 

0.242 

0.136 

-0.042 

0.795 

0.391 

0.136 

0.795 

m.I.e. 

24 

-0.074 

0.650 

-0.116 

0.475 

0.032 

0.846 

0.657 

0.475 

0.846 

m.I.e. 

25 

0.147 

0.364 

-0.253 

0.119 

-0.053 

0.746 

0.410 

0.119 

0.746 

m.I.e. 

26 

-0.042 

0.795 

0.021 

0.897 

-0.053 

0.746 

0.813 

0.746 

0.897 

m.m.m 

27 

-0.074 

0.650 

0.147 

0.364 

0.105 

0.516 

0.510 

0.364 

0.650 

m.m.m 

26 

-0.042 

0.795 

-0.263 

0.105 

-0.084 

0.604 

0.501 

0.105 

0.795 

m.I.e. 

26 

-0,074 

0.650 

-0.105 

0.516 

-0.105 

0.516 

0.561 

0.516 

0.650 

m.I.e. 

30 

-0.095 

0.559 

0.032 

0.846 

0.032 

0.846 

0.750 

0.559 

0.846 

m.I.e. 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

-0.08804 

0.45767 

-0.04211 

0.38100 

-0.01714 

0.55823 

0.46563 

0.24566 

0.69414 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

0.11499 

0.22407 

0.19334 

0.29522 

0.12607 

0.23159 

0.15545 

0.18088 

0.18533 
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Replication  Number 


Table  10.  Kendall  -Tau  Test  - 100  Engines 


KENDALL-TAU  TEST  OF  INDEPENDENCE  FOR  THE  SET  OF  100  ENGINES 


tau12 

pvalue 

tau13 

pvalue 

tau23 

pvalue 

mean 

min 

max 

best  method 

T 1 

0.029 

0.668 

-0.136 

0.045 

-0.072 

0.292 

0.335 

0.045 

0.668 

m.I.e. 

'2 

-0.093 

0.173 

-0.181 

0.008 

0.003 

0.967 

0.382 

0.008 

0.967 

m.I.e. 

■  3 

-0.181 

0.008 

0.003 

0.967 

-0.025 

0.708 

0.561 

0.008 

0.967 

m.I.e. 

4 

-0.144 

0.034 

-0.792 

0.428 

-0.064 

0.347 

0.270 

0.034 

0.428 

m.m.m 

5 

-0.089 

0.190 

0.088 

0.192 

-0.036 

0.592 

0.325 

0.190 

0.592 

m.m.m 

:!  e 

-0,138 

0.042 

-0.166 

0.015 

-0.098 

0.150 

0.069 

0.015 

0,150 

m.m.m 

..  7 

-0.028 

0.681 

-0.059 

0.381 

-0.064 

0.344 

0.469 

0.344 

0.681 

m.I.e. 

-0.036 

0.592 

-0.053 

0.435 

-0.145 

0.033 

0.353 

0.033 

0,592 

m.m.m 

ill 

-0.025 

0.708 

-0.009 

0.896 

-0.109 

0.108 

0.570 

0.108 

0.896 

m.I.e. 

fe  ill 

-0.026 

0.703 

-0.110 

0.104 

-0.019 

0.775 

0.527 

0.104 

0.775 

m.m.m 

1  -11 

-0.002 

0.976 

0.042 

0.536 

0.026 

0.703 

0.738 

0.536 

0.976 

m.I.e. 

1 

-0.132 

0.052 

0.064 

0.347 

-0.017 

0.798 

0.399 

0.052 

0.798 

m.I.e. 

§ 

-0.181 

0.008 

0.003 

0.967 

-0.025 

0.708 

0.561 

0.008 

0.967 

m.I.e. 

-0.081 

0.231 

-0.069 

0.308 

-0.018 

0.789 

0.443 

0.231 

0.789 

m.I.e. 

.a  i5 

-0.051 

0.449 

0.095 

0.162 

-0.093 

0.173 

0.261 

0.162 

0.449 

m.m.m 

S’  ”16 

-0.250 

0.000 

-0.060 

0.378 

-0.039 

0.568 

0.315 

0.000 

0.568 

m.I.e. 

^  17 

-0.070 

0.303 

0.011 

0.868 

-0.158 

0.020 

0,397 

0.020 

0.868 

m.I.e. 

•IS 

-0.059 

0.385 

0.016 

0.816 

-0.004 

0.953 

0.718 

0.385 

0.953 

m.m.m 

n 

-0.074 

0.273 

-0.053 

0.435 

-0.134 

0.048 

0.252 

0.048 

0.435 

m.I.e. 

-0.055 

0.415 

-0.022 

0.748 

-0.023 

0.739 

0.634 

0.415 

0,748 

m.I.e. 

.'21 

-0.058 

0.394 

-0.120 

0.077 

-0.068 

0.314 

0.262 

0.077 

0.394 

m.m.m 

■'22 

-0.104 

0.126 

0.122 

0,072 

-0.198 

0.004 

0.067 

0.004 

0.126 

m.I.e. 

-0.046 

0.501 

-0.161 

0.018 

-0.135 

0.046 

0.188 

0.018 

0.501 

m.I.e. 

:24 

-0.093 

0.173 

-0.118 

0.082 

-0.058 

0.391 

0.215 

0.082 

0.391 

m.I.e. 

0,022 

0.748 

0.071 

0.297 

-0.039 

0.568 

0.538 

0.297 

0.748 

m.I.e. 

-0.031 

0.651 

0.041 

0.544 

-0.210 

0.002 

0.399 

0.002 

0.651 

m.I.e. 

27 

0.006 

0,929 

-0.051 

0.457 

-0.002 

0.976 

0.787 

0.457 

0.976 

m.I.e. 

.28 

-0.029 

0.668 

-0.005 

0.938 

0.017 

0.803 

0.803 

0.668 

0.938 

m.m.m 

2$ 

-0.160 

0.019 

-0.022 

0.743 

-0.020 

0.770 

0.511 

0.019 

0.770 

m.I.e. 

■3d 

-0.189 

0.005 

0.002 

0.981 

-0.154 

0.024 

0.337 

0.005 

0.981 

m.m.m 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

mean 

-0.07884 

0.37010 

-0.05429 

0.44144 

-0.06601 

0.45688 

0.42281 

0.14569 

0.69136 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

stdev 

0.06784 

0.30264 

0.16079 

0.33413 

0.06379 

0.33785 

0.19490 

0.18499 

0.24739 
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V.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY 

The  analytical  approach  taken  to  develop  this  thesis  and  topics  covered  has  been  the 
following: 

1 .  Develop  a  program  in  S-Plus  to  generate  the  failure  times  for  an  engine  with  3 
subsystems,  where  the  distribution  of  the  first  failure  times  of  these  subsystems  is 
Weibull  and  the  engine  is  not  restored  to  new  condition  after  each  engine  failure. 

2.  Estimate  the  parameters  for  the  Weibull  distribution  that  fits  the  generated  engine 
failure  times  with  the  Maximum  Likelihood  Estimator  and  a  Modified  Method  of 
the  Moments  using  commercial  spread-sheet  software. 

3 .  Measure  the  accuracy  of  these  two  methods  with  the  Kolmogorov-Smimov  Test. 

4.  Apply  the  same  analysis  to  a  given  set  of  real  engine  failure  times. 

5.  Compare  the  accuracy  of  the  two  methods,  and  discuss  reasons  for  their 
difierences. 


B.  CONCLUSIONS 

The  two  methods  in  general  give  good  results  when  they  are  applied  to  fit  a  Weibull 
distribution  to  the  failure  times  generated  in  the  simulation.  The  resulting  fit  for  these  failure 
times  was  better  with  the  Maximum  Likelihood  Estimator  more  often  than  Avith  the  Modified 
Moment  Method.  For  this  reason,  and  the  fact  that  the  MLE  method  has  superior  theoretical 
properties  in  the  presence  of  independence,  the  author  recommends  the  MLE  method  to 
estimate  the  parameters  for  scenarios  similar  to  those  simulated  in  this  Thesis. 

With  these  conclusions  it  is  important  to  point  out  the  following  considerations: 
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The  MMM  has  the  advantage  that  fewer  calculations  are  necessary  to  get  the  results. 
In  the  case  of  the  Macros  developed  in  this  study,  only  the  initial  values  and  two  columns  of 
calculations  are  necessary.  Four  columns  of  calculations  are  necessary  in  the  MLE  case. 

There  was  a  wide  range  of  values  for  the  parameters  of  the  fitted  distributions  across  all 
the  replications.  This  is  displayed  in  Table  11. 


Table  11.  Parameters  Range 


MLE 

0.9611  <=  shape  <=  1.5238 
223.28  <=  scale  <=330.53 

20  Engines 

MMM 

0.9857  <=  shape  <=1.4616 

236.97  <=  scale  <=329.43 

MLE 

100  Engines 

MMM 

1. 1002  <=  shape  <=1,3422 

1.1737  <=  shape  <=1.4485 

255.64  <=  scale  <=301.25 

254.07  <=  scale  <=305.76 

The  dependence  obtained  in  the  simulation  is  not  as  significant  as  was  expected. 


C.  RECOMMENDATIONS 

These  two  methods  should  be  studied  when  greater  dependence  exists  among  the  data 
and  when  the  dependency  is  introduced  in  a  different  manner  than  used  in  this  thesis. 
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APPENDIX  A.  PROGRAM  TO  GENERATE  FAILURE  TIMES 


function(n) 

{ 

#  generate 

#  This  function  generates  times  randomly  in  a  manner  that  resembles  the  way  in  which 

#  engines  failures  occur  in  the  real  world. 

totalreplic  <-  NULL 
for(count  in  l:n)  { 

#  Step  1 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  first  failure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faultsl. 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  second  feilure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faults2. 

#  Three  Weibull  random  numbers  are  generated,  each  one  represents  the  third  failure  time 

#  that  each  one  of  the  components  could  have.  These  values  constitute  the  vector  faults3. 

xll<-rweibull(l,  1.2,  900) 
xl2  <-  rweibull(l,  1.3,  800) 
xl3  <-rweibull(l,  1.4, 600) 
faultsl  <«  c(xll,  xl2,  xl3) 
x21<-rweibull(l,  1.2,  900) 
x22<-rweibuU(l,  1.3, 800) 
x23  <-rweibull(l,  1.4, 600) 
faults2  <-  c(x21,  x22,  x23) 
x31<-rweibull(l,  1.2, 900) 
x32<-rweibuU(l,  1.3,800) 
x33  <-  rweibull(l,  1.4, 600) 
faults3  <-c(x31,x32,x33) 

failures  <-  matrix(c(faultsl,  faults2,  faults3),  ncol  =  3) 

#  Step  2 

#  The  minimum  value  of  the  vector  faultsl  is  selected,  this  value  is  the  time  of  the  first 

#  failure  of  the  engine.  The  program  creates  a  new  vector  newfail2  that  is  the  difference 

#  between  the  vector  faultsl  and  the  minimum  value  selected,  it  is  clear  that  this  new  vector 

#  has  a  component  that  is  0,  this  component  is  replaced  in  the  vector  newfail2  for  the 

#  element  of  the  vector  faults2  that  match  with  the  component  which  had  the  first  failure. 

#  Choose  the  first  failure  time,  that  is  the  minimum 

#  time  from  the  vector  faultsl 

fl  <-  min(faultsl) 

#  the  next  step  is  to  choose  the  second  failure  time,  to  make  this  first 

#  the  above  value  min(faultsl)  must  be  substituted  for  a  new  time 

#  generated  in  the  vector  faults2  and  precisely  in  the  same  component. 

#  After  this  the  new  minimum  is  selected 

newfail2  <-  (faultsl  -  fl) 
ii  <-  match(0,  newfail2) 
newfail2[ii]  <-  faults2[ii] 

#  Step  3 

#  At  this  point  the  engine  has  two  components  that  have  accumulated  working  hours  and 

#  one  that  has  been  just  repaired.  The  vector  newfail2  has  two  elements  that  have  the 
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#  remainder  of  the  theoretical  first  failure  for  these  components  and  one  element  that  is  the 

#  theoretical  second  failure  for  the  other  component.  So,  again,  the  minimum  element  of 

#  this  vector  is  selected  and  it  is  the  second  failure  time  of  the  engine. 


f2  <-  min(newfail2) 

#Step4 

#  The  program  creates  a  new  vector  called  newfaiB  that  is  the  difference  between  the  vector 

#  newfail2  and  this  minimum  value  selected,  it  is  clear  that  this  vector  has  a  component  that 
#isO. 

#  Choose  the  third  failure  time,  here  there  are  two  possible 

#  cases,  that  the  second  failure  is  in  the  same  component  that 

#  the  first  one  or  not,  these  two  condition  are  programmed  in  this  order 

# 

newfail3  <-  (newfail2  -  f2) 
iii  <-  match(0,  newfaiO) 
if(ii  =  iii){ 

newfail3[iii]  <-faults3[iii] 
f3  <-  min(newfail3) 

} 

else  if(ii  !=  iii)  { 
newfaiB  [iii]  <- faults2[iii] 

D  <“  min(newfail3) 

} 

enginefails  <-  c(fl,  f2,  f3) 

#cumenginefails  <-  cumsum(enginefails) 

#engfails  <-  matrix(c(failures,  enginefails,  cumenginefails),  byrow  = 

#  T,  ncol  =  3) 

totalrephc  <-  ibind(totalreplic,  (enginefails)) 

} 

totalreplic 

} 
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Row  Names 

VI 

V2 

V3 

1 

88.980852932 

7.374101636 

166.612433434 

2 

57.770945753 

452.164077047 

671.533220087 

3 

348.363054836 

13.725217444 

33.719418711 

4 

382.257870284 

168.795667771 

119.352260671 

5 

57.145359689 

117.795729216 

228.544288266 

6 

251.488400993 

417.099715694 

206.169150245 

7 

703.336492627 

16.866142294 

168.320579106 

8 

85.174113550 

960.378668893 

561 .250692879 

9 

126.379607253 

166.962891751 

410.881188239 

10 

9.871149363 

30.987527462 

1171.995737937 

11 

184.016586935 

55.183363300 

170.715797827 

12 

19.085994032 

206.177991742 

124.656084035 

13 

171.708369596 

7.411273439 

458.366290751 

14 

957.877021060 

334.941605556 

252.553811069 

15 

88.483061291 

129.532846338 

55.052226186 

16 

97.662694658 

64.764644477 

521 .223680534 

17 

213.546020387 

8.955316669 

314.364020817 

18 

420.340591055 

182.036946622 

345.506361499 

19 

46.823610226 

411.857588415 

109.450083602 

20 

25.300667331 

361.149199434 

744.757337329 

Column  1 

Column  2 

Column  3 

Column  1 

1 

Column  2 

1 

Column  3 

0.177625785 

1 
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Row  Names 

VI 

V2 

V3 

1 

450.634670353 

■fcgimgkfiE 

1  763.602013650 

2 

395.264314931 

300.696009190 

3 

467.906036969 

301 .510351673 

4 

200.004796884 

261 .419709923 

234.019837022 

5 

62.187250449 

293.954533661 

6 

188.420555520 

53_0.441439163 

446.499347177 

7 

517.045829232 

212.650323103 

4.372290612 

8 

104.406621033 

91 .440692508 

286.231002289 

9 

185.580658025 

10 

278.805459607 

142.687073420 

63.785901768 

11- 

242.374076871 

140.003133652 

57.571537504 

12 

35.982496826 

135.685578893 

13 

115.651429115 

807.515319374 

56.217725189 

14 

318.738490247 

220.511514630 

96.958868093 

15 

389.324028552 

121.699393963 

16 

833.642731237 

587.010089214 

17 

647.283140422 

270.006356536 

600.593877204 

18 

135.277879814 

235.260576261 

294.610909130 

19 

167.565552866 

80.687152199 

20 

43.940840941 

348.339229735 

21 

22 

198.697751449 

23 

24 

42.753769287 

568.105938757 

25 

290.367131217 

46.018913125 

470.135356798 

26 

442.572648090 

27 

315.869668084 

112.068427109 

799.073849390 

28 

258.473619251 

740.012832046 

730.391096774 

29 

936.074568617 

525.344669063 

53.458836253 

30 

38.275293391 

109.033459770 

31 

200.399299969 

564.453719467 

32 

16.387049764 

689.920420726 

124.595174314 

33 

. 282.346963237 

302.451537677 

17.105846707 

34 

298.005804703 

792.704752916 

277.592370478 

35 

488.070794225 

600.835486551 

^ilckl.g!W1l 

36 

651 .631129729 

186.384979015 

302.942800237 

37 

488.356688854 

282.631614706 

12.094239822 

38 

251 .322489919 

639.047414057 

315.553538531 

39 

102.245119335 

121 .065856465 

248.600567386 

40 

120.292743494 

281.418670027 

6.512039507 

41 

329.862437764 

42 

275.823361976 

400.377074953 

205.279868273 

43 

89.822151642 

568.897212277 

126.945901458 

44 

__  413.319725633 

280.829982419 

74.621453179 

45 

172.347938766 

138.081515648 

169.690890193 

46 

295.648493372 

47 

427.883569557 

120.358845277 

86.140070778 

48 

106.453003164 

336.696295308 

392.945334315 

49 

439.747178352 

145.848786206 

50 

95.424425759 

158.394437081 

233.2148472481 
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Row  Names 

VI 

V2 

V3 

51 

78.318115166 

233.292392917 

150.556841342 

52 

703.176710554 

131.985760999 

53 

284.082094482 

700.398323845 

54.502213913 

54 

673.813459742 

628.345709323 

27.810511050 

55 

416.157442911 

61.701311275 

56 

502.727434228 

208.984463635 

7.000471997 

57 

395.249167889 

181.615779362 

246.169600960 

58 

532.055233005 

282.774380172 

59 

256.526079351 

265.073575146 

366.480188363 

60 

791.013013348 

483.679778572 

383.122947523 

61 

116.674678096 

203.429431816 

501.860086251 

62 

216.019262031 

50.783576381 

346.268316228 

63 

297.710549841 

64 

57.143580125 

444.863422850 

53.260846841 

65 

305.430041853 

64.522458661 

344.328881540 

66 

134.875969833 

153.848996778 

74.573546508 

67 

207.795329813 

197.474623424 

36.572283433 

68 

148.585617029 

153.485311966 

69 

332.356937878 

183.167123218 

2.771444009 

70 

33.292174799 

281.412899135 

521.055842681 

71 

329.237256924 

624.659923834 

410.145804201 

72 

134.999172015 

924.566774196 

959.801158180 

73 

310.283780658 

196.949468983 

74 

145.927518733 

98.091275956 

105.358948583 

75 

96.240971550 

448.577731504 

68.349212862 

76 

504.457185330 

3.995392340 

54.819812986 

77 

779.185987694 

63.103630524 

510.654006981 

78 

64.295309451 

227.558484_64_1 

681.717084473 

79 

80.090112060 

305.523448635 

80 

225.006666940 

354.405022103 

115.791651829 

81 

322.833702315 

15.741543588 

273.716914308 

82 

552.256221658 

238.880493230 

103.236641646 

83 

267.774447132 

384.884639091 

84 

241.821803683 

27.660947307 

35.559153707 

85 

360.121368725 

54.400309193 

74.757012406 

86 

109.375105806 

169.046509653 

306.549361255 

87 

211.877188712 

47.413167726 

88 

328.701690509 

233.519929404 

411.204384691 

89 

271.788277670 

132.777355320 

90 

54.715071977 

52.466467632 

253.216046834 

91 

157.314280817 

49.880943701 

746.100150202 

92 

93.608619218 

110.699690363 

573.390687287 

93 

218.283430910 

434.282741182 

109.443120744 

94 

1106.554817220 

107.752314459 

95 

179.703059410 

96.191152623 

63.123472845 

96 

220.753490182 

358.476353401 

438.979104011 

97 

97.987549536 

57.556717482 

152.972644520 

98 

126.992796611 

173.350825362 

1004.747325381 

99 

656.574251587 

250.152057881 

31 2.989740534 

100 

189.330584986 

142.057479837 

63.629475350 

1  Cofumn  1  Column  2  Column  3  1 

11 

1 

0.012345873 

1 

-0.151985402 

0.027275574 

1 
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APPENDIX  B.  MACRO  20  ENGINES 


MACRO  FOR  20  ENGINES 
thesis  Macro 

Macro  recorded  by  FERNANDO  GUIMIL 


Sub  thesisO 
Sheets.Add 

Sheets("Sheet2”).Select 

Range("Al:C2").Select 

ActiveCelLFormulaRlCl  =  "MAXIMUM  LIKELIHOOD  ESTIMATOR" 

Range("AI:C2").Select 

Range("A2").Activate 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAligiunent  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xHorizontal 
End  With 
With  Selection 

.HorizontalAligtunent  =  xlCenterAcrossSelection 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xIHorizontal 
End  With 
Range("Dl").Select 

ActiveWindow.SmallScroll  ToRight:=10 
With  Toolbars(8) 

.Left  =  241 
.Top  =  228 
End  With 

Active  Window.LargeScroll  ToRight:=-l 
Range("Pl:Rl").Select 

ActiveCeU.FormulaRlCl  =  "MODIFIED  MOMENT  METHOD" 

Range("Pl:Rl").Select 

Range("QI").Activate 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xIHorizontal 
End  With 
With  Selection 

.HorizontalAhgnment  =  xlCenterAcrossSelection 
.VerticalAligiunent  =  xlBottom 
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.  WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 

RangeC’Nl").  Select 

ActiveWindow.LargeScrolI  ToRight:=-l 
WithTooIbars(8) 

.Left  =  436 
.Top  =  226 
End  With 
Range(”A2”).SeIect 

ActiveCell.FormulaRlCl  =  "shape  parameter" 

With  ActiveCell.Characters(Start:=l,  Length:=15).Font 
.Name  =  "Arial” 

.FontStyle  =  "Bold" 

.Size  =  10 

.Strikethrough  =  False 
.Superscript  =  False 
.Subscript  =  False 
.OutlineFont  =  False 
.Shadow  =  False 
.Underline  =  xlNone 
.Colorindex  =  xlAutomatic 
End  With 
Range("B3 ").  Select 

ActiveCell.FormulaRlCl  =  "initial  beta" 
Range("B3").Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Range("B5").Select 

ActiveCell.FormulaRlCl  =  "final  beta” 

Range("B5  ”).  Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  ==  xlHorizontal 
End  With 

Range("A7:C7").Select 
ActiveCell.FormulaRlCl  =  "scale  parameter" 
Range("A7:C7").Select 
Range("B7").Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
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.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("A7").SeIect 
ActiveCell.FonnulaRlCl  = "" 
Range("B7").Select 

ActiveCelLFomiulaRlCl  =  "scale  parameter" 

Range("B8").Select 

ActiveCelLFormulaRlCl  =  "landa" 

Range("B8").Select 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

ActiveWindow.SmallScroll  ToRight:=0 

Range("A2:C2").Select 

Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=l  1 
Range("Q2").Select 
ActiveSheetPaste 
Application.CutCopyMode  =  False 
Selection.Clear 

ActiveWindow.LargeScroll  ToRight;=-2 

Range("A2;C2").Select 

Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=9 
Range("P2").Select 
ActiveSheetPaste 
Range("04").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
Range("B3").Select 
Application.CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmaUScroll  ToRight:=ll 

Range("Q3").Select 

ActiveSheetPaste 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("B5").Select 
Apphcation.CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=12 

Range("Q5").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

ActiveSheetPaste 

Range("R7").Select 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("B5").Select 


Application.CutCopyMcxle  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=12 

Range("Q5").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight:=-2 
Range(”  A7 :  C8 ") .  Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=ll 
Range(”P7").Select 
ActiveSheet.Paste 
Range("N7”).Select 

ActiveWindow.LargeScroll  ToRight:=-2 

RangeC’Dll”).Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

RangeC'D7").Select 

ActiveSheet.Paste 

RangeC'BS").  Select 

ActiveWindow.SmallScroll  ToRight:=13 

ActiveWindow.LargeScroll  ToRight:=-2 

ActiveWindow.SmallScroll  ToRight:==14 

ActiveWindow.LargeScroU  ToRight:=-2 

RangeC’El").Select 

Application.  CutCopyMode  =  False 

ActiveCell.FormulaRlCl  =  "Number  of  observations" 

Range("El").Select 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAhgnment  =  xlCenterAcrossSelection 
.VerticalAlignment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 
Range("E2").Select 
ActiveCell.FormulaRlCl  -  "n" 

Range("E2").Select 

WithToolbars(8) 

.Left  =  536 
.Top  =  23 
End  With 

Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
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.VerticalAligiunent  =  xlBottom 
•WrapText  =  False 
.Orientation  =  xHorizontal 
End  With 

Sheets("Sheetl").Select 

Range("B2:B21").Select 

Selection.  Copy 

Sheets("Sheet2").Select 

Range("D10").Select 

ActiveSheetPaste 

Sheets("Sheetl").Select 

Range("C2:C21").Select 

Application.CutCopyMode  =  False 

SelectioaCopy 

Sheets("Sheet2").Select 

Range("E10").Select 

ActiveSheetPaste 

Sheets("Sheetl’').Select 

Range("D2:D21").Select 

Application.  CutCopyMode  =  False 

Selection.Copy 

Sheets("Sheet2'').Select 

Range("F10").Select 

Application.CutCopyMode  =  False 

Selection.Copy 

SheetsC'Sheetl  ").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheet2'').Select 

ActiveSheetPaste 

Range("D10:D29").Select 

AppUcation.CutCopyMode  =  False 

Selection.Copy 

Range("G10").Select 

ActiveSheetPaste 

Range("E10:E29").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

ActiveWindow.SmallScroll  Down:=16 

Range("G30").Select 

ActiveSheet.Paste 

Active  Window.LargeScroll  Down;=-l 
Range("F10:F29").Select 
Application.  CutCopyMode  =  False 
SelectioaCopy 

ActiveWindow.SmallScroll  Down;=35 

Range("G50").Select 

ActiveSheet.Paste 

Range("D42").Select 

ActiveWindow.LargeScroll  Down:=-2 

Range("G10:G69’’).Select 


Application.CutCopyMode  =  False 
Selection.  Copy 
Range(”H10").Select 
ActiveSheet.Paste 
Application.CutCopyMode  =  False 

Selection.Sort  Keyl:=Range("H10"),  Orderl:=xlAscending,  Header: 
xlGuess,  OrderCustom:=l,  ]VIatchCase:=False,  Orientation:=  _ 
xlTopToBottom 
Range('’B10").Select 
ActiveCell.FormulaRlCl  =  "theta" 

Range("B10").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAhgnment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Range("E3").Select 
ActiveCeU.FomiiilaRlCl  =  "60" 

Range("E3").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

ActiveWindow.ScrollRow=  1 
Range("Gl").Select 

ActiveCell.FormulaRlCl  =  'Tailures  times  Ti" 

Range("Gl").Select 
With  Selection 

.HorizontalAhgnment  =  xlCenterAcrossSelection 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  Trae 
.Orientation  =  xlHorizontal 
End  With 

SelectioaFonLBold  =  Tme 
Range("G9").Select 
ActiveCell.FormulaRlCl  =  "Ti" 

Range("G9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
SelectioaCopy 
Range("H9").Select 
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ActiveSheetPaste 
Range("G9").Select 
Application.  CutCopyMode  =  False 
Selection.  Clear 

ActiveWindow.SmallScroll  ToRight:=3 
Range("I9").Select 
ActiveCell.FormulaRlCl  =  "LnTi" 

Range(''I9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("I10").Select 

ActiveCell.FormulaRlCl  =  "=LN(RC[-1])" 
Range("I10").Select 

Selection.AutoFill  Destination:=Range("I10:I69"),  Type 
xlFiUDefault 
Range("I10:169").Select 
Range(’'H70").Select 
ActiveCeU.FormulaRlCl  =  "SUM" 
Range("H70").Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Selection.  Copy 
Range("I70").Select 
ActiveSheetPaste 
Range("H71").Select 
Application.CutCopyMode  =  False 
ActiveCelLFormulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 
Range("I71").Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-61]C:R[-21C)" 

Range("I72").Select 

ActiveWindow.LargeScroll  Down:=-3 

Range("J9").Select 

ActiveCeU.FormulaRlCl  =  "Ti^B4" 

Range("J9").Select 
Selection.Font.Bold  =  Trae 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
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ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroll  ToRight:=2 
Range('’J10”).Select 

ActiveCelLFormulaRlCl  =  ”-RC[-2]^R[-6]C[-8]’' 

Range("  J 1 0 ").  Select 

Selection.  AutoFill  Destination:=Range("J10:J69"),  Type:=  _ 
xlFillDefault 
Range("JiO:J69”).Select 

Selection.Replace  What:="^B*",  Replacement:="^B4",  LookAt:=xlPart, 
SearchOrder:=xlByRows,  MatchCase:=False 
Range(”K9”).Select 

ActiveWindow.SmallScroll  ToRight:=2 
ActiveCell.FormulaRlCl  =  'Ti^B6” 

Range("K9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("K10").Select 

ActiveCeU.FormulaRlCl  =  "=RC[-3]^R[4]C[-9]" 

RangeC'K  10").  Select 

Selection.  AutoFill  Destination:=Range("K10:K69"X  Type:=  _ 
xlFillDefault 

Range("K10:K69").Select 

Selection.Replace  What:="^B*",  Replacement:="^B6",  LookAt:=xlPart, 
SearchOrder:=xlByRows,  MatchCase:=False 
Range("L9").Select 

ActiveCelLFormulaRlCl  =  "Ln  TI*Ti'^B6" 

Range("L9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 
.VerticalAlignment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 
Range("L9").Select 

ActiveCelLFormulaRlCl  =  "Ln  TI  *  Ti^B6" 

With  ActiveCell.Characters(Start:=l,  Length:=13).Font 
.Name  = "  Arial" 

.FontStyle  =  "Bold" 

.Size=  10 

.Strikethrough  =  False 
.  Superscript  =  False 
.Subscript  =  False 
.OutlineFont  =  False 
.Shadow  =  False 
.Underline  =  xlNone 
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•Colorlndex  =  xlAutomatic 
End  With 

Range("L10").Select 

ActiveCeU.FonnulaRlCl  =  "=RC[-3]*RC[-1]" 
Range("L10").Select 

Selection.  AutoFill  Destination:=Range("L10:L69"),  Type 
xlFillDefeult 

Range("L10:L69").Select 
Range("M64").Select 
ActiveWindow-LargeScrolI  Down:=-3 
ActiveWindow.SmallScroll  Down:=12 
ActiveWindow.ScrollRow  =  28 
ActiveWindow.ScrollRow  =  1 
ActiveWindow.ScrollRow  =  61 
Range("K70").Select 
ActiveCelLFormulaRlCl  =  "SUM" 

Range("K70").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAligiunent  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Selection.Copy 
Range("L70").Select 
ActiveSheetPaste 
Range("K71").Select 
Application.CutCopyMode  =  False 
ActiveCeU.FormulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 
Range("L71").Select 

ActiveCelLFormulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 

Range("L72").Select 

With  Toolbars(8) 

.Left  =  57 
.Top  =  106 
End  With 

ActiveWindow.LargeScroll  Down:=-4 
Range("M9").Select 
ActiveCelLFormulaRlCl  =  "Ti'^" 

Range("M9").  Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAligrunent  =  xlCenter 
.  VerticalAligmnent  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("M10").Select 

ActiveCelLFormulaRlCl  =  "=RC[-5]''R[-4]C[4]" 
Range("M10").Select 
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Selection.  AutoFill  Destination:=Range("M10:M69''),  Type:=  _ 
xlFillDefault 

Range("M10:M69").Select 

SelectioaReplace  What:="^*",  Replacenient:="'XJ6",  LookAt:=xlPart, 
SearchOrder:=xlByRows,  MatchCase:=False 
ActiveWindow.SmallScroll  Down:=8 
Active  Window.  ScrollRow  =  33 
ActiveWindow.ScrollRow  =  43 
ActiveWindow.ScrollRow  =  48 
ActiveWindow.ScrollRow  =  51 
ActiveWindow.ScrollRow  =  57 
Range("M70").Select 
ActiveCell.FomiulaRlCl  =  "SUM" 

Range("M70").Select 
Selection.FontBold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("M71").Select 

ActiveCell.FomiulaRlCl  =  "=SUM(R{-61]C:R[-2]C)" 

Range("M72").Select 

ActiveCell.FomiulaRlCl  =  "mean" 

Range("M73").Select 
ActiveCell.FormulaRlCl  = "" 

Range("M72").Select 
Selection.Font.Bold  =  True 
With  Selection 

•HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
•WrapText  =  False 
•Orientation  =  xlHorizontal 
End  With 

Range("M73")^Select 
ActiveCeU^FormulaRlCl  =  "=R[-2]C/60" 

Range("M74")^Select 
ActiveWindow^ScroURow  =  1 
Active  Window^SmallScroU  ToRight;=2 
Range("N9")^Select 

ActiveCeU^FormulaRlCl  =  "(Ti'^-mean)''2" 

Range("N9")^Select 
With  Selection 

•HorizontalAlignment  =  xlCenterAcrossSelection 
•VerticalAlignment  =  xlBottom 
•WrapText  =  True 
•Orientation  =  xlHorizontal 
End  With 

SelectionFontBold  =  Trae 
Active  Window^SmallScroll  Down;=59 
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ActiveWindow-LargeScroll  Down:=4 
Range("N10").Select 

ActiveCeU.FormulaRlCl  =  "=(RC[-l]-R[63]C[-l]r2" 
Range("N10").Select 

Selection.  AutoFill  Destination~Range("N10;N69"),  Type:=  _ 
xlFillDefault 

Range("N10;N69").Select 

Selection.Replace  What:="-M*)^2",  Replacement:="-M73)^2",  LookAt 
:=xIPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("09").  Select 

ActiveCeU.FormulaRlCl  =  "(Ln  11)^2" 

Range("09").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("O10").Select 
ActiveCeU.FonnulaRlCl  =  "=RC[-6]^2" 

Range("O10").Select 

Selection.  AutoFiU  Destination:=Range("O10:O69"),  Type;=  _ 
xlFillDefault 

Range("O10;O69").Select 

Range("M70").Select 

Selection.  Copy 

Range("N70").Select 

ActiveSheetPaste 

Range("O70'').Select 

ActiveSheetPaste 

Range("N71").Select 

ApplicatioaCutCopyMode  =  False 

ActiveCeU.FonnulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 

Range("071").Select 

ActiveCelLFormulaRlCl  =  "=SUM(R[-61]C:R[-2]C)" 
Range("072").Select 
ActiveWindow.SmaUScroU  ToRight:=2 
ActiveWindow.LargeScroll  Down:=-3 
ActiveWindow.ScrollRow  =  1 
ActiveWindow.LargeScroU  ToRight:=-l 
WithToolbars(8) 

.Left  =  543 
.Top  =  37 
End  With 
Range("B4").Select 

ActiveWindow.SmaUScroU  ToRighl:=6 
ActiveWindow.SmaUScroU  Down:=60 
ActiveWindow.LargeScroU  Down:=4 
ActiveWindow.LargeScroU  ToRight;=-l 
Range("B4").  Select 
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ActiveCell.FormulaRlCl  =  _ 

"=((6/(PI0''2))*(R[67]C[13]-((R[67]C[7]^2)/60))*(l/59))^(-O.5)" 
Range("B4").Select 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=10 

Range("Q4").Select 

ActiveSheetPaste 

Application.  CutCopyMode  =  False 

ActiveCell.FonnulaRlCl  =  _ 

''=((6/(PI0^2))*(R[67]C[-2]-((R[67]C[-8]^2)/60))*(1/59))^(-0.5)" 

Range("Q5").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow-SmallScroll  ToRight:=-l 
RangeC’C6").Select 
ActiveCell.FormulaRlCl  =_ 
"=(60/RC[-l])+R[65]C[6]-(60*R[65]C[9]/R[65]C(8])" 
Range("C7").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
Range("B9").Select 

ActiveCellFormulaRlCl  =  "=(60/R[62]C[9])''(l/R[-3]C)" 

Range("B10").Select 

ActiveWindow.SmallScroll  ToRight:=-l 

Range("BH").Select 

ActiveCeU.FormulaRlCl  =  "=1/R[-2]C" 

Range("B12").Select 
ActiveWindow.SmallScroll  ToRight.-!  1 
Range("P6").Select 

ActiveCeU.FormulaRlCl  =  "=(R[65]C(-3])/((R[65]C[-2])^(0.5))" 
Range("Q9").Select 

ActiveWindow.LargeScroll  ToRight:=-2 
ActiveWindow.SmallScroll  ToRight:=9 
Range("Q10").Select 
ActiveCell.FormulaRlCl  =  "theta" 

Range("Q10").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
Range("Q9").Select 

ActiveCell.FormulaRlCl  =  "=1/(((1/60)*R[62]C[-4])''(1/R[-3]C))" 

Range("Qll").Select 

ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 

Range("Q12").Select 
End  With 


End  Sub 
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APPENDIX  C.  MACRO  100  ENGINES 


MACRO  TO  100  ENGINES 
thesis  Macro 

Macro  recorded  by  FERNANDO  GUIMIL 


Sub  thesisO 
Sheets.Add 

Sheets("Sheet2").Select 

Range("Al:C2").Select 

ActiveCeU.FormulaRlCl  =  "MAXIMUM  LIKELIHOOD  ESTIMATOR" 

Range("Al:C2").Select 

Range("A2").Activate 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAligrunent  =  xlCenter 
.  VerticalAligiunent  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 

Range("Dl").Select 

ActiveWindow.SmallScroll  ToRight:=10 
WithToolbars(8) 

.Left  =  241 
.Top  =  228 
End  With 

ActiveWindow.LargeScroll  ToRight:=-l 
Range("Pl:Rl").Select 

ActiveCelLFormulaRlCl  =  "MODIFIED  MOMENT  METHOD" 

Range("Pl:Rl").Select 

RangeC'Ql ").  Activate 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAligrunent  =  xlCenter 
.VerticalAligiunent  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAligiunent  =  xlCenterAcrossSelection 
.VerticalAlignment  =  xlBottom 
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.  WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 
RangeC'Nr’).Select 

ActiveWindow.LargeScroll  ToRight:=-l 
WithToolbars(8) 

.Left  =  436 
.Top  =  226 
End  With 
Range(''A2”).Select 

ActiveCell.FomiiilaRlCl  =  "shape  parameter" 

With  ActiveCell.Characters(Start:=l,  Length:=15).Font 
.Name  =  "Arial" 

.FontStyle="Bold" 

.Size=  10 

.Strikethrough  =  False 
.Superscript  =  False 
.Subscript  =  False 
.OutlineFont  =  False 
.Shadow  =  False 
.Underline  =  xlNone 
.Colorindex  =  xlAutomatic 
End  With 
Range("B3").Select 

ActiveCell.FormulaRlCl  =  "initial  beta" 
Range("B3").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Range("B5").Select 

ActiveCell.FormulaRlCl  =  "final  beta" 
Range("B5").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xIBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("A7:C7").Select 
ActiveCell.FormulaRlCl  =  "scale  parameter" 
Range("A7:C7").Select 
Range("B7").Activate 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
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.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("A7").SeIect 
ActiveCelLFormulaRlCl  = "" 
Range("B7").Select 

ActiveCell.FonnulaRlCl  =  "scale  parameter" 

Range("B8").Select 

ActiveCell.FonnulaRlCl  =  "landa" 

Range("B8").Select 

SelectLon.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAUgnment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

ActiveWindow.SmallScroll  ToRight:=0 

Range("A2:C2").Select 

Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=ll 
Range("Q2").Select 
ActiveSheetPaste 
ApplicatioaCutCopyMode  =  False 
Selection.Clear 

ActiveWindow.LargeScroll  ToRight;=-2 

Range("A2:C2").Select 

Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=9 
Range("P2").Select 
ActiveSheet.Paste 
Range("04").Select 

ActiveWindow.LargeScroll  ToRight:=-l 
Range("B3").Select 
Application.CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight~ll 

Range("Q3").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("B5").Select 
Application.CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=12 

Range("Q5").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

ActiveSheet.Paste 

Range("R7").Select 

ActiveWindow.LargeScroll  ToRight:=-2 
Range("B5").Select 


Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=12 

Range("Q5”).Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll  ToRight:=-2 
RangeC’A7:C8").Select 
Application.  CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  ToRight:=l  1 
Range("P7") .  Select 
ActiveSheet.Paste 
RangeC’N7”).Select 

ActiveWindow.LargeScroll  ToRight:=“2 

RangeCDir’).Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

RangeC’D7").Select 

ActiveSheet.Paste 

Range(”B8").Select 

ActiveWindow.SmallScroll  ToRight:=13 

ActiveWindow.LargeScroll  ToRight:=-2 

ActiveWindow.SmallScroll  ToRight:=14 

Active  Window.LargeScroll  ToRight:=-2 

Range(*’Er‘).Select 

Application.  CutCopyMode  =  False 

ActiveCell.FormulaRlCl  =  ’’Number  of  observations” 

Range(”El”).Select 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 
.  VerticalAhgnment  =  xlBottom 
.WrapText  =  Tme 
.Orientation  =  xlHorizontal 
End  With 
Range(”E2”).Select 
ActiveCell.FormulaRlCl  =  ”n” 

Range(”E2”).Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
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WithToolbars(8) 

.Left  =  539 
.Top  =  30 
End  With 

Sheets("Sheetl").Select 

Range("B2;B101").Select 

Selection.  Copy 

Sheets("Sheet2").SeIect 

ActiveWindow.SmallScroll  Down:=l 

Range("D10").Select 

ActiveSheetPaste 

Sheets("Sheetl").Select 

Range("C2;C101").Select 

Application.  CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheet2").Select 

Range("E10").Select 

ActiveSheet.Paste 

Sheets(''Sheetl").Select 

Range("D2:D101").Select 

Application.CutCopyMode  =  False 

Selection.  Copy 

Sheets("Sheet2").Select 

Range("F10").Select 

ActiveSheetPaste 

ActiveWindow.LargeScroll  ToRight:=0 
Range("D10:D109").Select 
Application.CutCopyMode  =  False 
Selection.  Copy 
Range("G10").Select 
ActiveSheetPaste 
Range("E10:E109").Select 
Application.CutCopyMode  =  False 
Selection.  Copy 

ActiveWindow.SmallScroll  Down;=100 

Range("G110").Select 

ActiveSheetPaste 

ActiveWindow.LargeScroll  Down:=-6 
Range("F10;F109").Select 
Application.  CutCopyMode  =  False 
Selection.  Cojty 

ActiveWindow.SmallScroll  Down:=200 

Range("G210").Select 

ActiveSheetPaste 

ActiveWindow.LargeScroll  Down:=-ll 
Range("G10:G309").Select 
Apphcation.  CutCopyMode  =  False 
Selection.  Copy 
Range("H10'').Select 
ActiveSheetPaste 
Application.CutCopyMode  =  False 


Selection.  Sort  Keyl:=Range("H10"),  Orderl:=xlAscendmg,  Header: 
xlGuess,  OrderCustom:=l,  MatchCase;=False,  Orientation:=  _ 
xlTopToBottom 

ActiveWindow.SmallScroll  Down:=-l 

Range("H2:H6").Select 

ActiveCell.FomiuIaRlCl  =  "failures  times  Ti" 

Range("H2:H6").Select 

Range("H3 ").  Activate 

Selection.Font.Bold  =  True 

With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 
.VerticalAlignment  =  xlBottom 
.WrapText  =  Tme 
.Orientation  =  xlHorizontal 
End  With 
Range("H9").Select 
ActiveCeU.FormulaRlCl  =  "Ti" 

Range("H9").Select 
Selection.FonLBold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Range("19").Select 
ActiveCeU.FormulaRlCl  =  "Ln  Ti" 

Range("I9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAligrunent  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("I10").Select 

ActiveCeU-FormulaRlCl  =  "=LN(RC[-1])" 

Range("I10").Select 

Selection.  AutoFiU  Destination:=Range("I10;I309"),  Type:=  _ 
xlFillDefeult 
Range("110:I309").Select 
WithToolbars(8) 

.Left  =  79 
.Top  =  270 
End  With 
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Range("G3 1 5").Select 
ActiveWindow.SmallScroll  ToRight:=2 
Range("I305 ").  Select 
ActiveWindow.LargeScroll  Down;=-17 
Range("J9").  Select 

ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroll  ToRight:=2 
Range("J9").  Select 
ActiveCeU.FormulaRlCl  =  "Ti^B4" 

Range("J9'').Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Selection.Copy 
Range("K9").Select 
ActiveSheet.Paste 

Active  Window.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroll  ToRight:=4 
Application.  CutCopyMode  =  False 
ActiveCeU.FonnulaRlCl  =  "Ti''B6" 

Range("J10").Select 

ActiveCeU.FormulaRlCl  =  "=(RC[-2]^R[-6]C[-8])" 
Range("J10").Select 

Selection.  AutoFiU  Destination:=Range("J10;J309"),  Type:=  _ 
xlFillDefault 
Range("J10:J309").Select 
Range("J300").Select 
ActiveWindow.LargeScroll  Down:=-16 
ActiveWindow.SmallScroll  Down;=-9 
Range("J10:J309").Select 

Selection.Replace  What:="^B*)",  Replacement:="'^B4)",  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("K10").Select 

ActiveCellFormulaRlCl  =  "=(RC[-3]'^R[-4]C[-9])" 
Range("K10").Select 

Selection.  AutoFill  Destination:=Range("K10:K309"),  Type;=  _ 
xlFillDefault 

Range(’’K10;K309").Select 

SelectioaReplace  What:="^B*)'',  Rq)lacement;="^B6)'',  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("L8:L9").Select 
ActiveCeU.FormulaRlCl  =  "Ln  Ti  *  Ti^B6" 
Range("L8:L9").Select 
Range("L9").  Activate 
With  Selection 

.HorizontalAlignment  =  xlCenterAcrossSelection 
.VerticalAhgnment  =  xlBottom 
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•WrapText  =  Trae 
■Orientation  =  xlHorizontal 
End  With 
Range("L9").Select 
Selection.DeIete  Shifl:=xlUp 
Range("L8").Select 
Selection.  Copy 
Range(”L9").Select 
ActiveSheet-Paste 
Range("L8").Select 
Application.  CutCopyMode  =  False 
Selection.  Clear 
Range(''L9").Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAhgnment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 

Range("L10").Select 

ActiveCellFormuIaRlCl  =  "=RC[-3]*RCI-1]" 
Range("L10").Select 

Selection.AutoFill  Destination:=Range("L10:L309"),  Type: 
xlFillDefault 

Range("L10:L309").Select 
ActiveWindow.LargeScroIl  Down:=-l 
Range("M285").Select 
ActiveWindow.LargeScroll  Down:=-16 
ActiveWindow.SmallScroU  ToRight:=4 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroU  ToRight:=8 
ActiveWindow.SmaUScroU  Down:=33 
Active  Window.LargeScroU  Down:=-2 
ActiveWindow.ScroURow  =  174 
Active  Window.  ScroURow  =  218 
ActiveWindow.ScroURow  =  252 
ActiveWindow.ScroURow  =  309 
Range("L310").Select 
ActiveCellFormuIaRlCl  =  "SUM" 

Range("L310").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAUgnment  =  xlCenter 
.VerticalAhgnment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("L311").Select 

ActiveCeU-FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 
Range("L312").Select 
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Range("M9'').Select 
ActiveCell-FormulaRlCl  =  "Ti'^" 

Range("M9").SeIect 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("M10").Select 

ActiveCell.FomiulaRlCl  =  "=(RC[-5]^R[-4]C[4])" 
Range("M10").Select 

Selection.  AutoFill  Destination:=Range("M10:M309"),  Type:=_ 
xlFillDefault 

Range("M10:M309").Select 

Selection.Replace  What:="'H3*)",  Replacement:="XJ6)",  LookAt: 

xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
Range("N9").Select 

ActiveWindow.SmallScroll  Down:=107 
ActiveWindow.ScrollRow  =  162 
AcliveWindow.ScrollRow  =  214 
ActiveWindow.ScrollRow  =  243 
ActiveWindow.ScrollRow  =  261 
ActiveWindow.ScrollRow  =  282 
ActiveWindow.ScrollRow  =  297 
Range("M310").Select 
ActiveCell.FomiulaRlCl  =  "SUM" 

Range("M3 10").  Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAhgnment  =  xIBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("M311").Setect 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 

Range("M312").Select 

ActiveCell.FormulaRlCl  =  "mean" 

Range("M312").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 
.VerticalAlignment  =  xIBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("M313").Select 
ActiveWindow.LargeScroll  Down;=-17 
Range("H3").Select 
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ActiveWindow.LargeScroll  ToRight:=-l 
Range("E3").Select 
ActiveCellFormulaRlCl  =  "300” 

Range(”E3").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.Vertical  Alignment  =  xIBottom 
,  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

ActiveWindow.SmallScroll  ToRight:=5 
ActiveWindow.ScrollRow=  144 
ActiveWindow,ScrolIRow=  178 
ActiveWindow.ScrollRow  =  219 
ActiveWindow.ScroURow  =  264 
ActiveWindow.ScrollRow  =  291 
ActiveWindow.ScrollRow  =  312 
ActiveWindow.ScrollRow  =  274 
ActiveWindow.LargeScroll  Down:=“15 
ActiveWindow.SmallScroll  Down:=37 
ActiveWindow.ScrollRow  =  283 
ActiveWindow.ScrollRow  =  312 
ActiveWindow.ScrollRow  =  303 
RangerM313”).Select 
ActiveCeU.FormulaRlCl  =  ”=R[-2]C/300” 
Range('’M314”).Select 
ActiveWindow.LargeScroll  Down:=-17 
Range(''N9”).  Select 

ActiveCell.FormulaRlCl  =  ”(Ti'X56“mean)^2” 
Range(”N9").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAhgnment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
With  Selection 

.HorizontalAlignment  =  xlCenter AcrossSelection 
.VerticalAlignment  =  xlBottom 
.WrapText  =  True 
.Orientation  =  xlHorizontal 
End  With 

RangeC‘N10”),Select 

ActiveCell.FormulaRlCl  =  ”=(RC[-2]-R[303]C[-l])” 
Range(”N10”).Select 

Selection.AutoFill  Destination:=Range("N10:N309”X  Type 
xlFLUDefault 

Range("N10:N309”).Select 
ActiveWindow.SmallScroll  ToRight:=l 
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RangerN10").Select 

ActiveCeU.FormiilaRlCl  =  "=(RC[-l]-RI303]C[-l]r2'‘ 

RangeC’N  10 ").  Select 

Selection.  AutoFill  Destination:=Range(”N10:N309"),  Type:=  _ 
xlFilDefault 

Range("N10:N309”).Select 

Selection.Replace  What:=*'-M*)^2'*,  Replacement  “”-M3 13)^2",  LookAt 
:=xlPart,  SearchOrder:=xlByRows,  MatchCase:=False 
ActiveWindow.SmallScroll  Down:=-2 
ActiveWindow.SmallScroll  ToRight:=3 
ActiveWindow.LargeScroll  ToRight:=-l 
RangeC'B4").Select 

ActiveCelLFonnulaRlCi  =  ''=((6/(PI0^2))*l)'’ 

Range("B5").Select 
ActiveWindow.SmallScroll  Down:=16 
ActiveWindow.ScrolIRow=  1 
ActiveWindow.SmallScroll  Down:=20 
Active  Window.  ScrolIRow  =  612 
ActiveWindow.ScrollRow  =  581 
ActiveWindow.  ScrolIRow  =  539 
ActiveWindow.ScrollRow  =  452 
ActiveWindow.LargeScroll  Down:=-9 
ActiveWindow.ScrollRow  =  339 
ActiveWindow.ScrollRow  =  309 
Range(”H310'’).Select 
ActiveCeU.FormulaRlCl  =  "SUM" 

Range("H310").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("H311").Select 

ActiveCelLFormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 
Range("H312").Select 
ActiveWindow.ScrollRow  =  1 
ActiveWindow.ScrollRow  =  465 
ActiveWindow.ScrollRow  =  287 
ActiveWindow.ScrollRow  =  322 
ActiveWindow.LargeScroll  ToRight:=l 
ActiveWindow.ScrollRow  =  203 
ActiveWindow.ScrollRow  =  269 
ActiveWindow.ScrollRow  =  3 
ActiveWindow.SmallScroll  ToRight:=4 
ActiveWindow.ScrollRow  =  124 
ActiveWindow.ScrollRow  =  223 
ActiveWindow.ScrollRow  =  317 
ActiveWindow.ScrollRow  =  274 
ActiveWindow.ScrollRow  =  276 
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ActiveWindow.SmallScroll  Down:=23 

RangeCH310”).Select 

Selection.  Copy 

Range("I310").Select 

ActiveSheet-Paste 

RangeC’I3ir’).Select 

Application.  CutCopyMode  =  False 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)” 

Range("I3 12").Select 

ActiveWindow.ScrollRow=  1 

ActiveWindow.SmallScroll  Down:=l 

ActiveWindow.SmallScroll  ToRight:=2 

Range(”09").Select 

ActiveCelLFormulaRlCl  =  ”(Ln  Ti)^2“ 

Range("09").Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAligmnent  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

RangeCOlO").  Select 
ActiveCell.FormulaRlCl  =  ’'=RC[-5]^2” 

Range("O10”).Select 

SelectioaAutoFill  Destination:=Range("010:0309"),  Type:= 
xlFillDefault 

Range("010:0309").Select 
Range("03 10").Select 
ActiveCell.FormulaRlCl  =  "SUM" 

Range("O310").Select 
Selection.Font.Bold  =  Tme 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("0311").Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 
Range("03 12").Select 
ActiveWindow.SmallScroll  ToRight:=-7 
ActiveWindow.LargeScroll  Down:=-15 
ActiveWindow.ScroilRow=  1 
ActiveWindow.SmallScroll  ToRight:=l 
Range("B4").Select 
ActiveCell.FormulaRlCl  =_ 

"=((6/(P10'^2))*(R[307]C[13]-R[307]C[7]/300)/299r(-0.5)" 

Range("B4").Select 

ActiveWindow.SmallScroll  ToRight:=6 
ActiveWindow.SmallScroll  Down:=-l 
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ActiveCell-FormulaRlCl  =_ 

"=((6/(PI0^2))*(R{307]C[13]-(RI307]C[7]^2)/300)/299)^(-0.5)" 
Range("B4").  Select 
ActiveWindow.ScrollRow  =  419 
ActiveWindow.ScrollRow  =  300 
ActiveWindow.SmallSCToll  ToRight:=7 
ActiveWindow.ScrollRow  =  1 
ActiveWmdow.LargeSaoll  ToRight:=-l 
ActiveCell.FormiilaRlCl  =_ 

"=((6/(PIO'^2))*(R[307]C[13]-(R[307]C[7]'^2)/300)/299r(-0.5)" 
Range("B4").Select 
ActiveCelLFormulaRlCl  =_ 

"=((6/(PI0^2))*(R[307]C[13]-(R[307]C[7]^2)/300)/299r(-0.5)" 
Range("I10").  Select 
Active  Window.SmallScroU  ToRight:=2 
Range("O10").Select 
Active  Window.LargeScroll  ToRight:=-l 
Range("B4").Select 
ActiveCelLFormulaRlCl  =  _ 

"=((6/(PI0^2))*(R[307]C[13]-(R[307]C[7]^2)/300)/299r(-0.5)" 
Range("B4").  Select 

ActiveWindow.SmallScroll  ToRight;=2 
Range("J10").Select 

ActiveWindow-LargeScroll  ToRight;=-l 
ActiveWindow.SmallScroll  ToRight:=2 
ActiveWindow.LargeScroll  ToRight:=-l 
Range("B4“).Select 

ActiveWindow.SmallScroll  ToRight:=2 
Range("Jll").Select 

ActiveWindow.LargeScroll  ToRight;=-l 
Range("B4").Select 
ActiveCelLFormulaRlCl  =  _ 

"=((6/(PI0^2))*(R[307]C[13]-(RI307]CI7]^2)/300)/299r(-0.5)" 
Range("O10").Select 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroll  ToRight:=7 
ActiveCeU.FormulaRlCl  =  "=RC[-ti]'^2" 

Range("O10").Select 

Selection.AutoFill  Destination:=Range("010:0309"),  Type:=  _ 
xlFillDefault 

Range("010:0309").Select 
ActiveWindow.LargeScroll  Down:=-16 
ActiveWindow.SmallScroll  Down:=-12 
ActiveWindow.LargeScroll  ToRight:=-l 
Range("B4").Select 

ActiveWindow.SmallScroll  ToRight:=3 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScroll  ToRight:=13 
Range("Q4").Select 

ActiveCelLFormulaRlCl  =  "=RC[-15]" 

Range("Q5”).Select 
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ActiveWindow.LargeScroll  ToRight:=-2 
WithToolbars(8) 

.Left  =  398 
.Top  =139 
End  With 

ActiveWindow.LargeScroll  ToRight:=l 
ActiveWindow.SmallScroll  ToRight:=-9 
Active  Window.SmallScroll  Down:=3 
Range(''C6").SeIect 
ActiveCell.FormulaRlCl  =_ 

"=(300/RC[-1])+R[305]C(6]-(300*R[305]C[9]/R[305]C[8])" 
Range("K310").Select 
ActiveCell.FormulaRlCl  =  "SUM" 

Range("K310").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAlignment  =  xBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("K311").Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 

Range("K312").Select 

ActiveWmdow.LargeScroll  Down:=-17 

ActiveWmdow.LargeScroll  ToRight:=-l 

Range("C6").Select 

ActiveCell.FormulaRlCl  =_ 

"=(300/RC[-1])+R[305]C[6]-(300*R[305]CI9]/R[305]C[8])" 

Range("C6").Select 

ActiveWindow.SmallScroll  ToRight:=2 
Active  Window.  ScrollRow  =  432 
ActiveWindow.ScrollRow  =  311 
Active  Window.  ScrollRow  =  1 
ActiveWindow.SmallScroll  ToRight:=-2 
Range("B10").Select 
ActiveCell.FormulaRlCl  =  "theta" 

Range("B10").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.VerticalAlignment  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 
Range("B9").Select 

ActiveCell.FormulaRlCl  =  "=(300/R[302]C[9]r(l/R[-3]C)" 

Range("B9").Select 

ActiveCeU.FormulaRlCl  =  "" 

Range("C9").Select 
ActiveCell.FormulaRlCl  = "" 
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Range("B9").Select 

ActiveCeU.FonnulaRlCl  =  "=(300/R[302]C[9])^(1/R[-3]C)" 

Range("Bll").Select 

ActiveCell.FormulaRlCl  =  "=1/R[-2]C" 

Range("Bll").Select 
ActiveWindow.SmallScrolI  ToRight:=9 
Range("Q6").Select 
WithToolbars(8) 

.Left  =  58 
•Top  =  25 
End  With 
Range("P6").Select 
ActiveCell.FormulaRlCl  =  "((1/299)" 

Range("S319").Select 
ActiveWindow.ScrollRow  =  302 
Range("N3 10").Select 
ActiveCell.FonnulaRlCl  =  "SUM" 

Range("N310").Select 
Selection.Font.Bold  =  True 
With  Selection 

.HorizontalAligiunent  =  xlCenter 
.  VerticalAlignment  =  xlBottom 
.  WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("N311").Select 

ActiveCell.FormulaRlCl  =  "=SUM(R[-301]C:R[-2]C)" 
Range("N312").Select 
ActiveWindow.LargeScroll  Down:=-17 
Range("P6").Select 

ActiveCell.FonnulaRlCl  =  "(M311)/((N311)^(0.5))" 
Range("P6").Select 

ActiveWindow.LargeScroll  ToRight:=0 
ActiveCell.FormulaRlCl  =  _ 

"=R[305]C[-3])/(R[305]C[-2])'^(0.5)" 

Range("Q9").Select 

ActiveCeU.FormulaRlCl  =  "=1/(((1/300)*RI302]C[-4])''(1/R[-3]C))" 
Range("Q10").Select 
ActiveWindow.LargeScroll  ToRight:=-l 
ActiveWindow.SmallScrolI  ToRight:=13 
ActiveCelLFormulaRlCl  =  "theta" 

Range("Q10").Select 
Selection.FonLBold  =  True 
With  Selection 

.HorizontalAlignment  =  xlCenter 
.  VerticalAligrunent  =  xlBottom 
.WrapText  =  False 
.Orientation  =  xlHorizontal 
End  With 

Range("Qll").Select 
ActiveCeU.FoimulaRlCl  =  "=1/R[-2]C" 
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Range("Q12").Select 
ActiveWindow.LargeScroll  ToRight:=-2 
End  Sub 
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APPENDIX  D.  PARAMETER  CALCULATIONS  FOR  THE  20  AND  100  ENGINES 


Row  Names  V1  V2  V3 


1 

88.980852932 

7.374101636 

166.612433434 

2 

57.770945753 

452.164077047 

671.533220087 

3 

348.363054836 

13.725217444 

33.719418711 

4 

382.257870284 

168.795667771 

119.352260671 

5 

57.145359689 

117.795729216 

228.544288266 

6 

251.488400993 

417.099715694 

206.169150245 

7 

703.336492627 

16.866142294 

168.320579106 

8 

85.174113550 

960.378668893 

561.250692879 

9 

126.379607253 

166.962891751 

410.881188239 

10 

9.871149363 

30.987527462 

1171.995737937 

11 

184.016586935 

55.183363300 

170.715797827 

12 

19.085994032 

206.177991742 

124.656084035 

13 

171.708369596 

7.411273439 

458.366290751 

14 

957.877021060 

334.941605556 

252.553811069 

15 

88.483061291 

129.532846338 

55.052226186 

16 

97.662694658 

64.764644477 

521.223680534 

17 

213.546020387 

8.955316669 

314.364020817 

18 

420.340591055 

182.036946622 

345.506361499 

19 

46.823610226 

411.857588415 

109.450083602 

20 

25.300667331 

361.149199434 

744.757337329 

1  Column  1  Column  2  Column  3  | 

Column  1 
Column  2 
Column  3 

1 

-0.099201099 

1 

-0.325682845 

0.177625785 

1 
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Row  Names 

VI 

V2 

V3 

1 

450.634670353 

134.620284266 

763.602013650 

2 

395.264314931 

300.696009190 

201.433889226 

3 

467.906036969 

301.510351673 

22.166715862 

4 

200.004796884 

261.419709923 

234.019837022 

5 

62.187250449 

6 

188.420555520 

530.441439163 

446.499347177 

7 

517.045829232 

212.650323103 

4.372290612 

8 

104.406621033 

91.440692508 

286.231002289 

9 

185.580658025 

189.658931128 

461.452317203 

10 

278.805459607 

142.687073420 

63.785901768 

11 

242.374076871 

140.003133652 

57.571537504 

12 

191.638384699 

35.982496826 

135.685578893 

13 

115.651429115 

807.515319374 

56.217725189 

14 

318.738490247 

220.511514630 

96.958868093 

15 

389.324028552 

16 

108.549643630 

833.642731237 

587.010089214 

17 

647.283140422 

270.006356536 

600.593877204 

18 

135.277879814 

235.260576261 

294.610909130 

19 

19.382052477 

80.687152199 

20 

620.417102151 

43.940840941 

348.339229735 

21 

549.652520297 

147.263663645 

402.329146256 

22 

160.858623029 

198.697751449 

131.936073250 

23 

216.407382425 

272.564254146 

24 

119.882883842 

42.753769287 

568.105938757 

25 

290.367131217 

46.018913125 

470.135356798 

26 

111.526180647 

laitMLifegiia 

442.572648090 

27 

112.068427109 

799.073849390 

28 

258.473619251 

740.012832046 

730.391096774 

29 

936.074568617 

525.344669063 

53.458836253 

30 

264.758222487 

109.033459770 

31 

200.399299969 

130.448227050 

564.453719467 

32 

16.387049764 

33 

302.451537677 

17.105846707 

34 

298.005804703 

792.704752916 

277.592370478 

35 

488.070794225 

600.835486551 

77.330402620 

36 

651.631129729 

186.384979015 

302.942800237 

37 

488.356688854 

282.631614706 

12.094239822 

38 

251.322489919 

39 

102.245119335 

40 

r-ettMiiiwiy-fj 

6.512039507 

41 

42.172683323 

42 

275.823361976 

400.377074953 

205.279868273 

43 

89.822151642 

568.897212277 

126.945901458 

44 

413.319725633 

280.829982419 

74.621453179 

45 

172.347938766 

138.081515648 

169.690890193 

46 

295.648493372 

49.777202258 

271.146965769 

47 

427.883569557 

120.358845277 

86.140070778 

48 

106.453003164 

336.696295308 

392.945334315 

49 

439.747178352 

145.848786206 

357.670884257 

50 

95.424425759 

158.394437081 

233.214847248 
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Row  Names 

VI 

V2 

V3 

51 

78.318115166 

233.292392917 

150.556841342 

52 

400.322169035 

703.176710554 

131.985760999 

53 

284.082094482 

lO 

00 

CO 

CM 

CO 

CO 

O) 

CO 

o 

o 

54.502213913 

54 

673.813459742 

628.345709323 

27.810511050 

55 

416.157442911 

61.701311275 

56 

502.727434228 

208.984463635 

7.000471997 

57 

395.249167889 

181.615779362 

246.169600960 

58 

532.055233005 

282.774380172 

51.530633142 

59 

256.526079351 

265.073575146 

366.480188363 

60 

791.013013348 

483.679778572 

383.122947523 

61 

116.674678096 

203.429431816 

501.860086251 

62 

216.019262031 

50.783576381 

346.268316228 

63 

652.055312871 

297.710549841 

64  • 

57.143580125 

444.863422850 

53.260846841 

65 

305.430041853 

64.522458661 

344.328881540 

66 

134.875969833 

153.848996778 

74.573546508 

67 

207.795329813 

36.572283433 

68 

55.079317773 

148.585617029 

153.485311966 

69 

332.356937878 

183.167123218 

2.771444009 

70 

33.292174799 

281.412899135 

521.055842681 

71 

329.237256924 

624.659923834 

410.145804201 

72 

134.999172015 

924.566774196 

959.801158180 

73 

310.283780658 

196.949468983 

61.420855422 

74 

145.927518733 

105.358948583 

75 

96.240971550 

448.577731504 

68.349212862 

76 

504.457185330 

3.995392340 

54.819812986 

77 

779.185987694 

63.103630524 

510.654006981 

78 

64.295309451 

227.558484641 

681.717084473 

79 

132.919226467 

80.090112060 

305.523448635 

80 

225.006666940 

354.405022103 

115.791651829 

81 

322.833702315 

15.741543588 

273.716914308 

82 

552.256221658 

238.880493230 

103.236641646 

83 

234.905588128 

267.774447132 

384.884639091 

84 

241.821803683 

27.660947307 

35.559153707 

85 

360.121368725 

54.400309193 

74.757012406 

86 

109.375105806 

169.046509653 

306.549361255 

87 

211.877188712 

47.413167726 

88 

328.701690509 

233.519929404 

89 

271.788277670 

132.777355320 

359.728088427 

90 

54.715071977 

52.466467632 

253.216046834 

91 

157.314280817 

746.100150202 

92 

93.608619218 

110.699690363 

573.390687287 

93 

218.283430910 

434.282741182 

109.443120744 

94 

1106.554817220 

178.424617471 

107.752314459 

95 

179.703059410 

96.191152623 

63.123472845 

96 

220.753490182 

fmmfm 

438.979104011 

97 

97.987549536 

57.556717482 

152.972644520 

98 

126.992796611 

173.350825362 

1004.747325381 

99 

656.574251587 

250.152057881 

312.989740534 

100 

189.330584986 

142.057479837 

63.629475350 

Column  1 

Column  2 

Column  3 

Column  1 

1 

Column  2 

0.012345873 

1 

Column  3 

<0.151985402 

0.027275574 

1 
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■ 

A  I  B...J _ C _ 

D 

E 

1 _ F  1  Q 

H 

\ 

J 

K 

L 

M 

N 

0 

P  \  Q  [  R 

1 

MAXIMUM  UKEUHOOD 
ESTIMATOR 

Number 

of 

observati 

ons 

MODIFIED  MOMENT  METHOD 

i 

shape  parameter 

n 

failures 
times  Ti 

shape  parameter 

3 

jinrtial  beta 

300 

[initial  beta 

4 

1.26241 

1.26241 

5 

final  beta 

final  beta 

6 

1.261028 

-1.9E-05 

17.34948 

1.288232 

m 

scale  parameter 

scale  parameter 

8 

landa 

landa 

9 

0.003418 

Ti 

LiiTi 

ri*B4 

Ti^BS 

LnTi* 

Ti*B6 

Ti*Q6 

(Ti*Q6- 

mean)*2 

(Ln  TD*2 

0.003394 

m 

theta 

450.6347 

134.6203 

763.602 

450.6347 

2.771444 

1.019368 

3.621402 

3.616305 

3.686347 

3.717989 

2290781 

1.039112 

theta 

EQ 

292.5569 

395.2643 

300.696 

201.4339 

395.2643 

3.995392 

1.385142 

5.746651 

5.735662 

7.944705 

5.955909 

2284012 

1.918618 

294.6595 

EQ 

467.906 

301.5104 

22.16672 

467.906 

4.372291 

1.475287 

6.439285 

6.426172 

9.480448 

6.689317 

2281796 

2.176472 

EO 

200.0048 

261.4197 

234.0198 

200.0048 

6.51204 

1.873653 

10.64743 

10.6199 

19.898 

11.17522 

2268263 

3.510574 

UL 

62.18725 

351.5663 

293.9545 

62.18725 

7.000472 

1.945978 

11.66534 

11.63401 

22.63953 

12.26648 

2264978 

3.786829 

EQ 

188.4206 

530.4414 

446.4993 

188.4206 

12.09424 

2.492729 

23.2626 

23.18261 

57,78798 

24.80917 

2227382 

6.213699 

EQ 

517.0458 

212.6503 

4.372291 

517.0458 

15.74154 

2.756303 

32.44627 

32.32293 

89.0918 

34.83971 

2197542 

7.597208 

Ea 

104.4066 

91.44069 

286.231 

104.4066 

16.38705 

2.796491 

34.13487 

34.00322 

95.08971 

36.69092 

2192057 

7.820364 

m 

185.5807 

189.6589 

461.4523 

185.5807 

17.10585 

2.83942 

36.03582 

35.89471 

101.9202 

38.77717 

2185884 

8.062308 

m 

278.8055 

142.6871 

63.7859 

278.8055 

19.14368 

2.951973 

41.53768 

41.36859 

122.1189 

44.82767 

2168030 

8.714144 

20 

242.3741 

140.0031 

57.57154 

242.3741 

19.38205 

2.964348 

42.19167 

42.0192 

124.5595 

45.54801 

2165909 

8.787356 

21 

191.6384 

35.9825 

135.6856 

191.6384 

22.16672 

3.098592 

49.98356 

49.77001 

1542169 

54.14712 

2140672 

9.601272 

22 

115.6514 

807.5153 

56.21773 

115.6514 

27.66095 

3.320022 

66.104 

65.80144 

218.4622 

72.02099 

2088689 

11.02254 

m 

318.7385 

220.5115 

96.95887 

318.7385 

27.81051 

3.325414 

66.55554 

66.25042 

220.3101 

72.52305 

2087238 

11.05838 

389.324 

121.6994 

391.9678 

389.324 

33.29217 

3.505322 

83.52575 

83.12216 

291.37 

91.43863 

2032940 

12.28728 

108.5496 

833.6427 

587.0101 

108.5496 

35.55915 

3.571198 

90.76888 

90.32207 

322.558 

99.53712 

2009912 

12.75345 

E 

6472831 

270.0064 

600.5939 

6472831 

35.9825 

3.583033 

92.13521 

91.66017 

328.493 

101 .0663 

2005578 

12.83812 

m 

1352779 

2352606 

294.6109 

135.2779 

36.57228 

3.599291 

94.04576 

93.57919 

336.8187 

1032054 

1999524 

12.95489 

28 

167.5656 

19.38205 

80.68715 

167.5656 

3827529 

3.644805 

99.60763 

99.10723 

3612265 

109.4375 

1981938 

13.2846 

29 

620.4171 

43.94084 

348.3392 

620.4171 

42.17268 

3.741773 

112.5787 

111.9981 

419.0714 

123.9987 

1941151 

14.00086 

30 

549.6525 

1472637 

402.3291 

549.6525 

42.75377 

3.755457 

114.5404 

113.9476 

427.9253 

126.204 

1935011 

14.10346 

31 

160.8586 

198.6978 

131.9361 

160.8586 

43.94084 

3.782844 

118.5697 

117.9516 

446.1925 

130.7361 

1922423 

14.30991 

32 

269.7486 

216.4074 

272.5643 

269.7486 

46.01891 

3.829052 

125.6921 

125.0288 

478.7419 

138.7547 

1900251 

14.66164 

E3 

119.8829 

42.75377 

568.1059 

119.8829 

47.41317 

3.8589 

130.5185 

129.8244 

500.9794 

144.1938 

1885285 

14.89111 

EH 

290.3671 

46.01891 

470.1354 

290.3671 

49.7772 

3.907557 

138.787 

138.0396 

539.3977 

153.5213 

1859758 

15.269 

35 

111.5262 

59.04016 

442.5726 

111.5262 

49.88094 

3.909639 

139.1522 

138.4025 

541.1038 

153.9337 

1858633 

15.28528 

36 

315.8697 

112.0684 

799.0738 

315.8697 

50.78358 

3.927573 

142.3386 

141.5682 

556.0193 

157.5314 

1848837 

15.42563 

37 

258.4736 

740.0128 

730.3911 

258.4736 

61.53063 

3.942176 

144.987 

144.1994 

568.4593 

160.523 

1840710 

15.54076 

38 

936.0746 

525.3447 

53.45884 

936.0746 

52.46647 

3.960174 

148.3189 

147.5095 

584.1633 

164.2883 

1830507 

15.68298 

39 

38.27529 

264.7582 

109.0335 

38.27529 

53.26085 

3.975201 

151.1594 

150.3314 

597.5976 

167.4997 

1821828 

15.80223 

200.3993 

130.4482 

564.4537 

200.3993 

53.45884 

3.978912 

151.8691 

151.0365 

600.9608 

168.3022 

1819662 

15.83174 

ED 

16.38705 

689.9204 

124.5952 

16.38705 

54.40031 

3.99637 

155.2534 

154.3984 

617.033 

172.1302 

1809349 

15.97097 

Ea 

282.347 

302.4515 

17.10585 

282.347 

54.42758 

3.996871 

155.3516 

154.496 

617.5005 

172.2413 

1809050 

15.97498 

m 

298.0058 

792.7048 

277.5924 

298.0058 

54.50221 

3.998241 

155.6206 

154.7632 

618.7806 

172.5457 

1808232 

15.98593 

ED 

488.0708 

600.8355 

77.3304 

488.0708 

54.71507 

4.002139 

156.3882 

155.5258 

622.4358 

173.4143 

1805896 

16.01712 

o 

651.6311 

186.385 

302.9428 

651.6311 

54.81981 

4.004052 

156.7663 

155.9013 

6242369 

173.842 

1804747 

16.03243 

46 

488.3567 

282.6316 

12.09424 

488.3567 

55.07932 

4.008774 

157.7037 

156.8325 

628.7062 

174.9029 

1801898 

16.07027 

47 

251.3225 

639.0474 

315.5535 

251.3225 

56.21773 

4.029232 

161.8296 

160.9311 

648.4288 

179.5736 

1789380 

16.23471 

48 

102.2451 

121.0659 

248.6006 

1022451 

57.14358 

4.045567 

1652014 

164.2805 

664.6077 

183.3925 

1779178 

16.36661 

49 

1202927 

281.4187 

6.51204 

1202927 

57.55672 

4.052771 

166.7106 

165.7796 

671.8669 

185.1023 

1774619 

16.42495 

SO 

394.6289 

329.8624 

42.17268 

394.6289 

57.57154 

4.053028 

166.7648 

165.8335 

672.1277 

185.1637 

1774456 

16.42704 

Ea 

275.8234 

400.3771 

2052799 

275.8234 

59.04016 

4.078218 

172.1531 

171.1857 

698.1326 

1912708 

1758223 

16.63186 

52 

89.82215 

568.8972 

126.9459 

89.82215 

61.42086 

4.117749 

180.9624 

179.9356 

740.9299 

201.2637 

1731822 

16.95586 

53 

413.3197 

280.83 

74.62145 

413.3197 

61.70131 

4.122305 

182.0061 

180.9723 

746.0232 

202.4484 

1728705 

16.9934 

54 

172.3479 

138.0815 

169.6909 

172.3479 

62.18725 

4.13015 

183.8176 

182.7715 

754.8737 

204.5047 

1723302 

17.05814 

55 

295.6485 

49.7772 

271.147 

295.6485 

63.10363 

4.144778 

1872436 

186.1743 

771.6513 

208.3951 

1713103 

17.17919 

427.8836 

120.3588 

86.14007 

427.8836 

63.12347 

4.145093 

187.318 

1862481 

772.0158 

208.4795 

1712882 

17.18179 

izi 

106.453 

336.6963 

392.9453 

106.453 

63.62948 

4.153077 

1892155 

188.1328 

781.33 

210.6348 

1707245 

17.24805 

ELI 

439.7472 

145.8488 

357.6709 

439.7472 

63.7859 

4.155532 

189.803 

188.7162 

7842163 

211.3021 

1705502 

17.26845 

ELI 

_ ! _ 

95.42443 

158.3944 

2332148 

95.42443 

64.29531 

4.163487 

191.7185 

190.6187 

793.6385 

213.4785 

1699822 

17.33462 

m 

78.31812 

2332924 

150,5568 

78.31812 

64.52246 

4.167013 

192.574 

191.4683 

797.8511 

214.4506 

1697286 

17.364 

lil 

400.3222 

703.1767 

131.9858 

400.3222 

68.34921 

4.22463 

207.103 

205.8976 

869.841 

230.9735 

1654509 

17.8475 

m 

284.0821 

700.3983 

54.50221 

284.0821 

74.57355 

4.311786 

231.1907 

229.8173 

990.923 

258.4184 

1584659 

18.5915 

m 

673.8135 

628.3457 

27.81051 

673.8135 

74.62145 

4.312428 

231.3782 

230.0035 

991.8735 

258.6323 

1564120 

18.59704 

EH 

416.1574 

416.1574 

230.5305 

994.56471  2592377 1 

1582597 

18.61269 
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APPENDIX  E.  FUNCTION  ^^tocompare” 


function(v,  betamle,  thetamle,  betammm,  thetamnun) 

{  #  This  function  calculates  the  kolmogorov-smimov  goodness  of  the  fit, 

#  and  plots  of  the  empical  distribution  and  theoretical  distribution 

#  for  a  set  of  data  and  two  sets  of  parameters  of  weibuU  distribution. 
win.graphO 

a  <-  ks.gof(v,  dist  =  "weibull",  shape  =  betamle,  scale  =  thetamle) 
b  <-  ks.gof(v,  dist  =  "weibull",  shape  =  betammm,  scale  =  thetamnun) 
par(m£row  =  c(l,  2)) 

cdf.compare(v,  dist  =  "weibull",  shape  =  betamle,  scale  =  thetamle) 
cdf  compare(v,  dist  =  "weibull",  shape  =  betammm,  scale  =  thetamnun) 
list(a,  b) 

} 
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tocompare(a20. 1,0.973, 25 1.71, 0.9857, 253. 14) 

[[1]]: 

One-sample  Kolmogorov-Smimov  Test;  hypothesized  distribution  =  weibull 

data;  v 

ks  =  0.0546,  p-value  =  0.994 

alternative  hypothesis:  Tme  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 

[[2]]: 

One-sample  Kolmogorov-Smimov  Test;  hypothesized  distribution  =  weibull 
data:  v 

ks  =  0.0509,  p-value  =  0.9977 

Alternative  hypothesis:  Tme  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 
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Empirical  and  Hypothesized  weibull  CDFs  Empirical  and  Hypothesized  weibull  CDFs 


-solid  iiife  h  Ihc  pirical  6f. 

MLE 


3oiid  lioe  n-^tKsompiriald/. 

mmm; 


Figure  5.  Kolmogorov-Smimov  Test;  Hypothesized  Distribution  =  Weibull.  (a)  ks  = 
0.0546,  p-value  =  0.994,  (b)  ks  =  0.0509,  p-value  =  0.9977 


tocompare(al00.1, 1.261, 292.56, 1.2882, 294.66) 

[[1]]: 

One-sample  Kolmogorov-Smimov  Test;  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.03 19,  p-value  =  0.9207 

Alternative  hypothesis:  True  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 

[[2]]: 

One-sample  Kolmogorov-Smimov  Test;  hypothesized  distribution  =  weibull 

data:  v 

ks  =  0.0382,  p-value  =  0.7741 

Alternative  hypothesis:  Tme  cdf  does  not  equal  the  weibull  Distn.  for  at  least  one  sample  point. 
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Empirical  and  Hypothesized  weibull  CDFs  Empirical  and  Hypothesized  weibull  CDFs 


'Mlid  lt*K  h  IhB  pircal  dA 


solid  line  isiheempiriald/. 


MLE 


MMM 


Figure  6.  Kolmogorov-Smimov  Test;  Hypothesized  Distribution  =  Weibull.  (a)  ks  = 
0.0319,  p-value  =  0.9207,  (b)  ks  =  0.0382,  p-value  =  0.7741 
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